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Chapter 1 
Introduction 



In this review I would like to show the importance and the power of measurements 
of anisotropies in the cosmic microwave background (CMB). 

CMB anisotropies are so useful mainly because they are small: For a given model, 
they can be calculated within linear perturbation theory, to very good approxima- 
tion. They are influenced only little by the non-linear processes of galaxy formation. 
This allows us to compute them very precisely (to about 1%, which is high preci- 
sion for present cosmological standards). For given initial fluctuations, the result 
depends only on the cosmological parameters. If we can measure CMB anisotropies 
to a precision of, say 1%, this allows us therefore to determine cosmological parame- 
ters to about 1%. An unprecedented possibility! Consider that at present, after the 
work of two generations, e.g. the Hubble parameter is known only to about 25%, 
the baryon density is known to about 10% and the uncertainties in the dark matter 
density, the cosmological constant and the space curvature are even larger. 

This somewhat too optimistic conclusion has however three caveats which we 
want to mention before entering the subject of this review. 

1. Initial conditions: The result depends on the model for the initial fluctu- 
ations. The simplest inflationary scenarios which lead to adiabatic perturba- 
tions, contain in general three to four free parameters, like the ratio of tensor to 
scalar perturbations (r) and the spectral index of the scalar and tensor pertur- 
bations {ns and n^), so a few more parameters need to be fitted additionally 
to the data. 

More generic initial conditions allow for at least four additional isocurvature 
modes with arbitrary (anti-)correlations. The initial conditions are then given 
by a 5 X 5 positive semi-definite matrix, and, in principle, several spectral in- 
dices 01 . In most of this review we shall ignore this possibility and assume 
that initial perturbations are purely adiabatic. Even if isocurvature constribu- 
tions cannot be excluded, this most simple model is in good agreement with 
the present data. 

If the perturbations are generated by active sources like, e.g., topological de- 



fects, then the modehng is far more comphcated, and the analysis is too dif- 
ferent to be included in this review. 

2. Degeneracy: Even though we can measure over 1000 independent modes 
(Cf's) of the CMB anisotropy spectrum, there are certain combinations of the 
cosmological parameters that lead to degeneracies in the CMB spectrum. The 
result is, e.g., very sensitive to the sum r2matter + ^a, but not to the difference 
( "cosmic confusion" ) . 

3. Cosmic variance: Since the fluctuations are created by random processes, 
we can only calculate expectation values. Yet we have only one universe to 
take measurements ("cosmic variance"). For small-scale fluctuations we can 
in general assume that the expectation value over ensembles of universes is the 
same as a spatial average (a kind of ergodic hypothesis), but for large scales 
we can't escape large statistical errors. 



1.1 Friedmann-Lemaitre universes 

Priedmann-Lemaitre universes are homogeneous and isotropic solutions of Einstein's 
equations. The hyper-surfaces of constant time are homogeneous and isotropic, i.e., 
spaces of constant curvature with metric a^{r))jijdx^dx^ , where 7,^ is the metric of 
a space with constant curvature k. This metric can be expressed in the form 

-fijdx'dx^ = dr^ + x\r) {d^^ + siri^'&dip'^) (1.1) 

X^(r) = <j sinV , k = 1 (1.2) 
y sinh^ r , K, — —1, 

where we have rescaled a{r)) such that k — ±1 or 0. (With this normalization the 
scale factor a has the dimension of a length and 77 and r are dimensionless ior k 0.) 
The four-dimensional metric is then of the form 



Qn^dx'^dx" = —a^{vi)drf' + a^{ri)'yijdx^dxK (1.3) 

Here 77 is called the conformal time. 

Einstein's equations reduce to ordinary differential equations for the function 
0(77) (with ■ = d/drf): 

-] +K = — -aV+-Aa' 1.4 
a J 6 6 
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where p = —Tq, p = (no sum!) and all other components of the energy momen- 
tum tensor have to vanish by the requirement of isotropy and homogeneity. A is the 
cosmological constant. 

Energy momentum "conservation" (which is also a consequence of ( |1.4| ) and ( p..5| ) 
due to the contracted Bianchi identity) reads 

P = -3(^)(P + P). (1.6) 

After these preliminaries (which we suppose to be known to the audience) let 
us answer the following question: Given an object with comoving diameter AQ at a 
redshift z{ri) = (ao/a) — 1. Under which angle ^{X,z) do we see this object today 
and how does this angle depend on and fi^? 

We define 




V=VO 

Aa2 



3(f)' 



(1.7) 

n=vo 



iff 



V=VO 



where the index o indicates the value of a given variable today. Friedmann's equation 
(|0|) then requires 

1 = fi^ + l^A + (1.8) 
Back to our problem: Without loss of generality we set r = at our position 



z=z„ ri=rij 

©---v/:;;;;;;;;_'_''']o' 



rio,z=0 - «A 



B 

Figure 1.1: The two ends of the object emit a flash simultaneously from A and B at zi 
which reaches us today. 



and thus r = ri = tjq — rji at the position of the flashes, A and B at redshift Zi. If A 
denotes the comoving arc length between A and B we have A = xi'^^i)^ = x{Vo~Vi)'^j 
i.e. 

^= , ^ V (1-9) 

xiVo - Vi) 



^or physical size a(?7)A — d 
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It remains to calculate {rjo — r]i){zi). 

Note that in the = we can still normalize the scale factor a as we want, 

and it is convenient to choose = 1, so that comoving scales today become physical 
scales. However, for k 7^ 0, we have already normalized a such that n = ±1 and 
X = sinr or sinhr. We have in principle no normalization constant left. 

From the Friedmann equation we have 

a^ = ^a''p+-Aa''-Ka\ (1.10) 
3 3 

We assume that p is a combination of "dust" (cold, non-relativistic matter) with 
Pd = and radiation with prad = VsPrad- 

From ( pTBI) we find that p^ad oc a""^ and pd oc a~^. Therefore, with Hq = (^) (rjo), 
we define 

SirG 



3 



aV = {a^Qr^d + Qaaal) (1.11) 



-Aa^ = H^Q^a^ (1.12) 
3 

-Ka^ = H^Q^al . (1.13) 
The Friedmann equation then implies 

da r, ( ^ a ^ ^ a"^ ^ \ ^ 

— = Hoal[n,^d + —^d + —^A + —^^] 1.14 

df] \ oo Oq / 

so that 

1 r dz 

7]0-7]l = — / -. (1.15) 

^0^0 io [n^^Z + ly + nd{z + 1)3 + l^A + ^.{z + 

Here we have introduced the cosmological redshift z + 1 = a^/a. (In principle we 
could of course also add other matter components hke, e.g. "quintessence" 0, which 
would lead to a somewhat different form of the integral ( |1.15| ), but for definiteness, 
we remain with dust, radiation and a cosmological constant.) 

In general, this integral has to be solved numerically. It determines the angle 
'(?(A, zi) under which an object with comoving size A at zi is seen. 

On the other hand, the angular diameter distance to an object of physical size d 
seen under angle 'd is given by rjQ — rji = ri = (of^) • ^^^^ measure 

the redshift and the comoving angular diameter distance of a certain class of objects 
comparing with Eq. (|1.15|) allows in principle to determine the parameters Qm, ^a, 

and Hq. 

We have = ^ H^aQ = for ^ 0. 



Observationally we know 10^^ < i^rad < 10^^ as well as 0.1 < ~ 1, \^a\ ~ 1 
and \VtA ~ 1. 
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n. = 0; n,= -0.8, -0.3, 0, 0.3, 0.1 
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Figure 1.2: The function x(% ~ ^i) a function of the redshift z for different values 
of the cosmological parameters (left, with f2A=0) and f^A (right, with 0^=0), namely 
-0.8 [dotted], -0.3 [short-dashed], [sohd], 0.3 [dot-dashed], 0.8 [long-dashed]. 



If we are interested in small redshifts, z\ ~ 10, we may safely neglect firad- In 
this region, Eq. ( |1.15|) is very sensitive to and provides an excellent mean to 
constrain the cosmological constant. 

At high redshift, z\ ~ 1000, neglecting radiation is no longer a good approxima- 
tion. 

We shall later need the opening angle of the horizon distance, 

^h{z,) = (1.16) 
Xivo - Vi) 

^OaoA, [^]^^^(^+l)4 + f]^(^ + l)3 + ^^ + ^^(^ + l)2]r 

(Clearly this integral diverges if flj-^d = = 0. This is exactly what happens during 
an inflationary period and leads there to the solution of the horizon problem.) 

The value of the radiation density is well known. For photons plus three sorts 
of massless neutrinos we have 

p,ad = 7.94 X 10-3^(To/2.737K)%/cm=^ . 

This gives 

n.^dh^ = 4.2 ■ 10-^(ro/2.737K)^ , (1.18) 

Ho = 100h^. (1.19) 
sMpc 

Neglecting Qra.d, for f^A = and small curvature, < [^2^1 ^ at high enough 



redshift, zi > 10, one has rjo — rji ~ 2^J\V^^l\/V^d = 2/(i^oOovm)- This yields 
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r!,= 0; = -0.8, -0.3, 0, 0.3, O.f 



n„ = 0; n.= -0.8, -0.3, 0, 0.3, o.f 
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Fi gure 1.3: '&Hizi) (in degrees) for different values of the cosmological parameters ri^j 



and the line styles are as in Fig. 1.2 



??(A, zi) ~ y/ri^HoaoX/2 = ^y/TUHoXphy s/{zi + 1), where Xphys = ctiA is the physical 
scale corresponding to comoving size A. 

1.2 Recombination and the cosmic microwave back- 
ground (CMB) 

During its expansion, the universe cools adiabatically. At early times, it is domi- 
nated by a thermal radiation background with p = C/a^ = gcso-^BT^^ and we find 
that T oc a~^. Here Qeff = + T/Sni;' is the effective number of degrees of freedom, 
bosons counting as 1 and fermions counting as 7/8 (see e.g. 0]). At temperatures 
below O.SMeV only neutrinos and photons are still relativistic leading to the density 
parameter given in Eq. ( |1.18| ). ( Neutrinos have a somewhat lower temperature than 
photons, = {^./llfl^T, since they have already dropped out of thermal equilib- 
rium at T ~ IMeV, before annihilation which therefore reheats the photons but 
not the neutrinos, see e.g. [|l^, p]].) 

The photons obey a Planck distribution. 



1 



(1.20) 



At a temperature of about T ~ 4000K ~ 0.4eV, the number density of pho- 
tons with energies above the hydrogen ionization energy drops below the baryon 



^We will use units with h — c = = 1 throughout this report. The Stefan-Boltzniann constant 
is then given by osb = ir^k^/ {GOfV^c^) — 7r^/60. 
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density of the universe, and the protons begin to (re-) combine to neutral hydrogen. 
(Hehum has aheady recombined earher.) Photons and baryons are tightly coupled 
before (re-) combination by non-relativistic Thomson scattering of electrons. During 
recombination the free electron density drops sharply and the mean free path of the 
photon grows larger than the Hubble scale. At the temperature Tjec ~ 3000K (corre- 
sponding to the redshift ^dec — HOO and the physical time t^ec = cioVdec — lO^years) 
photons become free and the universe becomes transparent. 

After recombination, the photon distribution evolves according to Liouville's 
equation (geodesic spray): 

p%f - Vp''^^ ^ LxJ = 0, (1.21) 

where i = 1,2, 3. Since the photons are massless, |pp = Yl^=iPiP^ = LiP' [uj = ap^). 
Isotropy of the distribution implies that / depends on only via |p| = uj, and so 

In a Friedmann universe (also if k 7^ 0!) we find for p'^p^ = — cu^ + = [exercise!] 

r»>. = -u^' f 4l • (1-23) 



Inserting this result into ( 1.21 ) leads to 



a\ df 



which is satisfied by an arbitrary function / = f{uja). Hence the distribution of 
free-streaming photons changes just by redshifting the momenta. Therefore, setting 
T (X even after recombination, the blackbody shape of the photon distribution 
remains unchanged. 

Note however that after recombination the photons are no longer in thermal 
equilibrium and the T in the Planck distribution is not a temperature in the thermo- 
dynamical sense but merely a parameter in the photon distribution function. 

The blackbody spectrum of these cosmic photons which are called the "cos- 
mic microwave background" (CMB) is extremely well verified observationally (see 
Fig. |T]^). The limits on deviations are often parameterized in terms of three param- 
eters: The chemical potential /x, the Compton y parameter (which quantifies a well 
defined change in the spectrum arising from interactions with a non-relativistic elec- 
tron gas at a different temperature, see e.g. lTl|]) and (describing a contamination 
by free-free emission). 

The present limits on these parameters are (at 95% CL, [0) 

|/i|<9■10~^ |^/|<1.2■10-^ iFffl < 1.9- 10"^ (1.25) 
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Figure 1.4: Spectrum of the cosmic background radiation. The graph on the left shows 
the measurements of the FIRAS experiment on COBE (the vertical bars), overlaid by a 
blackbody spectrum at a temperature of 2.73 K. The error bars are 20 times magnified! 
The image on the right shows a larger number of measurements. The FIRAS data is 
represented by the fat line around the peak of the spectrum (from Peebles Q]). 



The CMB Photons have not only a very thermal spectrum, but they are also 
distributed very isotropically, apart from a dipole which is (most probably) simply 
due to our motion relative to the surface of last scattering: 

An observer moving with velocity v relative to a source emitting a photon with 
proper momentum p = —un sees this photon redshifted with frequency 

cj' = (1 - nv) , (1.26) 

where 7 = is the relativistic 7-factor. For an isotropic emission of photons 

coming from all directions n this leads to a dipole anisotropy in first order in v. 
This dipole anisotropy, which is of the order of 

(-) 



~ 10" 

dipole 



has already been discovered in the 70ties ||T2|, |T^. Interpreting it as due to our 
motion with respect to the last scattering surface implies a velocity for the solar- 
system bary-center of v = 371 ± 0.5km/s at 68% CL (0)- 

The COBEQ DMR experiment (Differential Microwave Radiometer) has found 
fluctuations of 



10-^ (1.27) 



on all angular scales 6 > 7° [^. On smaller angular scales many experiments have 
found fluctuations (we shall describe the experimental results in more detail later), 
but all of them are ~ 10~^. 



^Cosmic Background Explorer, NASA satellite launched 1990. 
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As we shall see later, the CMB fluctuations on large scales provide a measure 
for the deviation of the geometry from the Friedmann-Lemaitre one. The geometry 
perturbations are thus small and we may calculate their effects by linear perturbation 
theory. On smaller scales, bT jT reflects the fluctuations in the energy density in 
the baryon/radiation plasma prior to recombination. Their amplitude is just about 
right to allow the formation of the presently observed non-linear structures (like 
galaxies, clusters, etc.) out of small initial fluctuations by gravitational instability. 

These findings strongly support the hypothesis which we assume here, namely 
that the large scale structure [i.e. galaxy distribution) observed in the universe 
formed by gravitational instability from relatively small (~ 10~^ — 10~^) initial fluc- 
tuations. As we shall see, such initial fluctuations leave an interesting "fingerprint" 
on the cosmic microwave background. 
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Chapter 2 

Perturbation Theory 



The tool for the analysis of CMB anisotropies is cosmological perturbation theory. 

We spend therefore some time on this subject, especially on the fundamental level. 

Once all the variables are defined, we will be rather brief in the derivation of the 
basic perturbation equations. First of all, because these derivations are in general 
not very illuminating and secondly because nowadays all of you can obtain them 
very easily by setting 

S-zxi/ = Qixu + ea^h^,, (2.1) 

{9nv being the unperturbed Priedmann metric) and asking Mathematic or Maple 
to calculate the Einstein Tensor using the condition — 0. We conventionally set 
(absorbing the "smallness" parameter e into h^i,) 

g^iu = g,iu + (i^h^,^, goQ^-a^, Qij^a^lij \h^^\<^l 

2.1 Gauge transformation, gauge invariance 

The first fundamental problem we want to discuss is the problem of 'choice of gauge' 
in cosmological perturbation theory: 

For linear perturbation theory to apply, the spacetime manifold 7W with metric 
g and the energy momentum tensor T of the real, observable universe must be in 
some sense close to a Friedmann universe, i.e., the manifold M. with a Robertson- 
Walker metric g and a homogeneous and isotropic energy momentum tensor T. 
It is an interesting, non-trivial unsolved problem how to construct g and T from 
the physical fields g and T in practice. There are two main difficulties: Spatial 
averaging procedures depend on the choice of a hyper-surface of constant time and 
do not commute with derivatives, so that averaged fields g and T will in general 
not satisfy Einstein's equations. Secondly, averaging is in practice impossible over 
super-horizon scales. 
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Even though we cannot give a constructive prescription, we now assume that 
there exists an averaging procedure which leads to a Friedmann universe with spa- 
tially averaged tensor fields Q, such that the deviations {Tfj,u — T^y)/ max{o/3}{|TQ^|} 
and {g^iu — g^u)/ Hiax{a/3}{^^^} are small, and g and T satisfy Friedmann's equations. 
Let us call such an averaging procedure 'admissible'. There may be many differ- 
ent admissible averaging procedures (e.g. over a different hyper-surface) leading to 
slightly different Friedmann backgrounds. But since \g — g\ is small of order e, the 
difference of the two Friedmann backgrounds must also be small of order e and we 
can regard it as part of the perturbation. 

We consider now a fixed admissible Friedmann background {g,T) as chosen. 
Since the theory is invariant under diffeomorphisms (coordinate transformations), 
the perturbations are not unique. For an arbitrary diffeomorphism and its pullback 
0*, the two metrics g and 4'*{g) describe the same geometry. Since we have chosen 
the background metric g we only allow diffeomorphisms which leave g invariant 
i.e. which deviate only in first order form the identity. Such an 'infinitesimal' 
isomorphism can be represented as the infinitesimal flow of a vector field X, (p = (p^ . 
Remember the definition of the fiow: For the integral curve 7a; (s) of X with starting 
point X, i.e., •y^is = 0) = x we have 4>f{x) = 7a;(s). In terms of the vector field X, 
to first order in e, its pullback is then of the form 

(p* = id + eLx 

{Lx denotes the Lie derivative in direction X). The transformation g <P*{g) is 
equivalent to ^ + ta?h ^ g + t{a?h + Lxg), i-e. under an 'infinitesimal coordinate 
transformation' the metric perturbation h transforms as 

h h + a~^Lxg . (2.3) 

In the context of cosmological perturbation theory, infinitesimal coordinate trans- 
formations are called 'gauge transformation'. The perturbation of a arbitrary tensor 
field Q = Q + eQ^^^ obeys the gauge transformation law 

Q(i) ^ g(i) + LxQ . (2.4) 

Since every vector field X generates a gauge transformation (p = 0^, we can 
conclude that only perturbations of tensor fields with LxQ = for all vector fields 
X, i.e., with vanishing (or constant) 'background contribution' are gauge invariant. 
This simple result is sometimes referred to as the 'Stewart Walker Lemma' 0. 

The gauge dependence of perturbations has caused many controversies in the lit- 
erature, since it is often difficult to extract the physical meaning of gauge dependent 
perturbations, especially on super-horizon scales. This has led to the development 
of gauge invariant perturbation theory which we are going to use throughout this 
review. The advantage of the gauge-invariant formalism is that the variables used 
have simple geometric and physical meanings and are not plagued by gauge modes. 
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Although the derivation requires somewhat more work, the final system of pertur- 
bation equations is usually simple and well suited for numerical treatment. We 
shall also see, that on sub-horizon scales, the gauge invariant matter perturbations 
variables approach the usual, gauge dependent ones. Since one of the gauge invari- 
ant geometrical perturbation variables corresponds to the Newtonian potential, the 
Newtonian limit can be performed easily. 

First we note that since all relativistic equations are covariant (i.e. can be written 
in the form Q = for some tensor field Q), it is always possible to express the 
corresponding perturbation equations in terms of gauge invariant variables 0, ||, . 

2.2 Gauge invariant perturbation variables 

Since the {i] = const} hyper-surfaces are homogeneous and isotropic, it is sensible 
to perform a harmonic analysis: A (spatial) tensor field Q on these hyper-surfaces 
can be decomposed into components with transform irreducibly under translations 
and rotations. All such components evolve independently. For a scalar quantity / 
in the case k = this is nothing else than its Fourier decomposition: 

/(x,r^) = I d\f{y)e''^-. (2.5) 

(The exponentials yk(x) = e*'''' are the unitary irreducible representations of the 
Euclidean translation group.) For n = 1 such a decomposition also exists, but the 
values k are discrete, k"^ = i{i + 2) and for k, = —1, they are bounded from below, 
k"^ > 1. Of course, the functions Ik are different for k ^ 0. 

They are always the complete orthogonal set of eigenfunctions of the Laplacian, 

AF(^) = -fcV(^). (2.6) 

In addition, a tensorial variable (at fixed position x) can be decomposed into 
irreducible components under the rotation group SO (3). 

For a vector field, this is its decomposition into a gradient and a rotation, 

Vi = Vi^ + Bi, (2.7) 

where 

B\, = 0, (2.8) 

where we used X\i to denote the three-dimensional covariant derivative of X. ip is 
the spin and B is the spin 1 component of V. 
For a symmetric tensor field we have 

H., = i?L7., + (v.V, - iA7.,) Ht + 1 (iff,;) + H^l^) + h\J\ (2.9) 
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where 

= hP' = H^' = 0. (2.10) 

(V) (T) 

Here and Ht are spin components, is a spin 1 component and H>j is a 
spin 2 component. 

We shall not need higher tensors (or spinors). As a basis for vector and tensor 
modes we use the vector and tensor type eigenfunctions to the Laplacian, 



= -evP (2.11] 

and 



Al^f = -k'rp , (2.12) 

where Yj^^ is a transverse vector, Y^.*-^-"-' = and Yji'^ is a symmetric transverse 
traceless tensor, Y^'^^^ = vjp^^ = 0. 

According to Eqs. (|2.?| ) and (|2.9|) we can construct scalar type vectors and tensors 
and vector type tensors. To this goal we define 

Yj-^^ = -k^^Y^f (2.13) 

^ fc-i^(f + 17.->^(^) (2.14) 

In the following we shall extensively use this decomposition and write down the 
perturbation equations for a given mode k. 

The decomposition of a vector field is then of the form 

Bi = BY}^^ + B^^^yI^I (2.16) 
The decomposition of a scalar field is given by (compare p.9| ) 

H,, = H,Y('h, + HtYIP + H^^^Y^P + H(^^Y}P, (2.17) 
where B, bP , Hi, Ht, hP and hP are functions of rj and k 



2.2.1 Metric perturbations 

Perturbations of the metric are of the form 

9f,u = gf,v + o^h^^. (2.18) 

We parameterize them as 

h^^dx^'dx" = -2Adrf - 2B,dridx' + 2Hijdx'dx\ (2.19) 
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and we decompose the perturbation variables Bi and Hij according to ( |2.16| ) and 

Let us consider the behavior of hf^i, under gauge transformations. We set the 
vector field defining the gauge transformation to 



X = TS^ + Udi. 

Using simple identities from differential geometry like L^{df) = d{L:^f) and 
{Lxl)ij = Xi\j + Xj\i, we obtain 



Ly^g 



^2.20) 



2\-T + T]dri' + 2[Li-TA drjdx 



+ I 2-T'jij + Liy + Ljii ) dx^dx^ 



:2.21) 



Comparing this with ( p.l9| ) and using (|2.3| ) we obtain the following behavior of 
our perturbation variables under gauge transformations (decomposing Li = LY-^ + 



A A + -T + T 
a 

B B-L-kT 

a 6 
Ht Hr-kL 



(2.22) 

(2.23) 
(2.24) 

(2.25) 

(2.26) 
(2.27) 
(2.28) 



Two scalar and one vector variable can be brought to disappear by gauge transfor- 
mations. 

One often chooses kL = Ht and T = B + L, so that the variables Ht and B 
vanish. In this gauge (longitudinal gauge), scalar perturbations of the metric are of 
the form {Ht = B = 0): 

h^s) = _2^ciry2 ^ 2^-fijdx'dx^. (2.29) 
\E' and $ are the so called Bardeen potentials. In general they are defined by 

^ = A--k-^(T-k-^& (2.30) 
a 

$ = Hr. + \ht - -k-^a (2.31) 



a 



with a = k ^Ht — B. A short calculation using Eqs. (|2.22|) to ( p.26| ) shows that 
they are gauge invariant. 
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For vector perturbations it is convenient to set kL^^'' = H^^^ so that H^^^ 
vanishes and we have 

h^Ux^'dx'' = 2a^^^YPd7]dx\ (2.32) 

We shall call this gauge the "vector gauge". In general a^^^ = k~^H^^^ — B^^^ is 
gauge invariantFl. 

Clearly there are no tensorial (spin 2) gauge transformation and hence Hl^ is 
gauge invariant. 

2.2.2 Perturbations of the energy momentum tensor 

Let T^^ = t;; + be the full energy momentum tensor. We define its energy density 
p and its energy flow 4- vector u as the time-like eigenvalue and eigenvector of Tj^: 

Tl^u" = -pu^, = -1. (2.33) 

We then define their perturbations by 

p = p(l + 5), u = u^dt + u'di. (2.34) 

M° is fixed by the normalization condition, 

M° = -(l-A). (2.35) 
d 

We further set 

u' = -v' = + v^^^Y^^^\ (2.36) 

a 

We define Pj^ = u'^Ui, + 6j^, the projection tensor onto the part of tangent space 
normal to u and set the stress tensor 

T^"' = P^P^T""^. (2.37) 

In the unperturbed case we have Tq = 0,rj = p6j. Including perturbations, to 
first order we still obtain 

r° = r° = = 0. (2.38) 
But rj contains in general perturbations. We set 

rj = p[{i + Ul) 6] + n;.] , with = o. (2.39) 

We decompose 11* as 

n} = n(^)Ff )^ + n(^)y/^)^ + (2.40) 

We shall not derive the gauge transformation properties in detail, but just state 
some results which can be obtained as an exercise (see also 0): 

the shear of the hyper-surfaces of constant time. 
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Of the variables defined above only the n^'^'^'^) are gauge invariant; they de- 
scribe the anisotropic stress tensor, 11^ = r/f — ^/^t^S'^- They are gauge in- 
variant due to the Stewart- Walker lemma, since n = 0. For perfect fluids 
W. = 0. 



A second gauge invariant variable is 



r = - ^S, (2.41) 

w 

where cj. = p/p is the adiabatic sound speed and w = p/p is the enthalpy. One 
can show that F is proportional to the divergence of the entropy flux of the 
perturbations. Adiabatic perturbations are characterized by F = 0. 

• Gauge invariant density and velocity perturbations can be found by combining 
5, V and f j-^"* with metric perturbations. 

We shall use 

V = v-\ht = v^^''''^^ (2.42) 

A/ 

Dg = 5 + 3(1 + w) [hl + = + 3(1 + w)<!> (2.43) 

D = + 3(1 +w)(-]^ (2.44) 

\a J k 

ViV) = W) = ^(vec) (2.45) 

n = _ = - (2.46) 
fi-l^W = a^^\ (2.47) 

Here v^^"''^^ and ff^*"^'' are the velocity (and density) perturbations in the 

longitudinal and vector gauge respectively and a(^) is the metric perturbation in 
vector gauge (see Eq. ( p.32| )). These variables can be interpreted nicely in terms of 



gradients of the energy density and the shear and vorticity of the velocity field |]14 
But we just want to show that on scales much smaller than the Hubble scale, 

ki] ^ 1, the metric perturbations are much smaller than S and v and we can thus 

"forget them" (which will be important when comparing experimental results with 

calculations in this formalism): 

The perturbations of the Einstein tensor are given by second derivatives of the 

metric perturbations. Einstein's equations yield the following order of magnitude 

estimate: 

'5T\ „ ,^ „ / 1 . 



O [ O {SttGTI = O \^-^h +'-h + k^hj (2.48) 

= 0{h + kr]h+{kr]fh) . (2.49) 
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For krj ^ 1 this gives 0{6,v) = O (^) ^ 0{h). On sub-horizon scales the differ- 
ence between 5, Dg and D is negligible as well as the difference between v 
and V or v(^), and l^(^). 

Later we shall also need other perturbation variables like the perturbation of the 
photon brightness (energy-integrated photon distribution function), but we shall 
introduce them as we get there and discuss some applications first. 

2.3 Basic perturbation equations 

As already announced, we do not derive Einstein's equations but just write down 
those which we shall need later: 



2.3.1 Constraint equations 



A-nGa^pD = (fc^ - 3ft;)$ (00) 
A'KGa?{p + p)V = k[{l)-^-^^ (Oi) 

%'KGa'{p + p)^l = \{2K-k^)a^^^ (Oi) 



(scalar) 
(vector) 



(2.50) 
(2.51) 



2.3.2 Dynamical equations 



-A;^($ + *) = STrGa^pn^^) (scalar) (2.52) 
= STrGa^pn^^) (vector) (2.53) 



fc(^W + 2(-)a(^) 
a 



if{^) + 2(^^^i7(^) + (2«; + A;2)//(^) = STrGa^pn^f (tensor) (2.54) 

There is a second dynamical scalar eqn., which is however complicated and not 
needed, since we may instead use one of the conservation eqns. below. Note that for 
perfect fluids, where IT*- = 0, we have $ = — ^, oc 1/a^ and B. obeys a damped 
wave equation. The damping term can be neglected on small scales (over short time 
periods) when rf'^ ~ 2k -|- A;^, and H^^ represents propagating gravitational waves. 
For vanishing curvature, these are just the sub-horizon scales, kr] ~ 1. For k < 0, 
waves oscillate with a somewhat smaller frequency, uj — \/2k-\- /c^, while for k > 
the frequency is somewhat larger. 



2.3.3 Conservation equations 



D, + 3 (c^ - w) (f) I?, + (1 + w)kV + ?>w (j) r 

wk 



+ 



l+w 



[r-i(i-|l)n] 



(scalar) (2.55) 
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W^A'-tY-' (vecto.) (2.56) 
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Chapter 3 

Simple applications 



We first discuss some simple applications which will be important for the CMB. We 
could of course also write ( |2.55| ) in terms of D, but we shall just work with the 
relation 



D = Dg + ?>{l + w)[-^+[-]k-W]. (3.1 



3.1 The pure dust fluid at /t: = 0, A = 

We assume the dust to have = = p = and 11 = F = 0. The equations ( p.55| ), 
( p.52| ) and (|2.50|) then reduce to 



Dg = —kV (energy conservation eqn.) (3.2) 



+ I - M/ = k'^ (gravitational acceleration eqn.) (3.3) 



$ = (3.4) 
-k^^ = A7cGa^p(^Dg + 3(^^ + k-^^^ (Poisson eqn.). (3.5) 

In a pure dust universe p oc =^ {a/aY oc a^^, which is solved by a oc rj"^. 
The Einstein equations then give immediately AnGpa? = 2/2(0/0)^ = G/rj"^. Setting 
ki] = X and ' = d/dx, the system ( p.2[ - p.5|) then becomes 

(3.6) 
(3.7) 

^ ^Dg + 3f^ + -F)) = (3.8) 



(3.9) 



\ " \ X 





-V 


V = 




X 









We use (p.8[) to eliminate \E' and (|3^) to eliminate Dg, leading to 



(18 + x') V" + + V -(^-^ + A^V = 
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The general solution of Eq. ( p.9|) is 

l^ = K)X+^ (3.10) 

with arbitrary constants Vq and Vi. Since the perturbations are supposed to be small 
initially, they cannot diverge for x — >■ 0, and we have therefore to choose Vi = (the 
growing mode). Another way to argue is as follows: If the mode Vi has to be small 
already at some early initial time r/in, it will be even much smaller at later times 
and may hence be neglected. The perturbation variables are then given by 

V = Vox (3.11) 

Dg = -15Vo-^VoX^ (3.12) 

^ = 3Vo. (3.13) 

The constancy of the gravitational potential \1/ in a matter dominated universe 
and the growth of the density perturbations like the scale factor a led Lifshitz to con- 



clude 1946 [|T5[ that pure gravitational instability cannot be the cause for structure 
formation: If we start from tiny thermal fluctuations of the order of 10~^^, they can 
only grow to about lO^'^'^ through this process during the matter dominated regime. 
Or, to put it differently, if we do not want to modify the process of structure for- 
mation, we need initial fluctuations of the order of at least 10~^. One possibility to 
create such fluctuations is due to quantum particle production in the classical grav- 
itational field during inflation. The rapid expansion of the universe during inflation 
quickly transforms microscopic scales at which quantum fluctuations are important 
into cosmological scales where these fluctuations are then "frozen in" as classical 
perturbations in the energy density and the geometry. 

We distinguish two regimes: 
i) super-horizon, a; -C 1 where we have 

Dg = -15Vo (3.14) 
^ = 3Vo (3.15) 
V = Vox (3.16) 

and ii) sub-horizon, x ^ 1 where the solution is dominated by the terms 

V = Vox (3.17) 

Dg = —Vox' (3.18) 

^ = 3V^) = const (3.19) 

6., 1, 



Note that for dust 



D = Dg + S"^ + -V = --Vox 
X 2 



2 



In the variable D the constant term has disappeared and we have D <^ on super 
horizon scales, x <^ 1. 
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3.2 The pure radiation fluid, a€ = 0, A = 



In this limit we set w — = 1/3 and H = 0. We conclude from p oc that a oc rj 
and $ = — ^, and the perturbation equations become (with the same notation as 
above) : 



v - 

The general solution of this system is 



2* + -Dg 



1 

4' 



12. 



3D„ + 12* + —V 

X 



Dr. 



Do 



cos 



f X \ ^/S ( X 
I — 1 - 9^ sin I — 



VV3 



X 



VV3 



sm 



f X \ ^\/3 f X 
1^1+ 2— cos I 



X 



V 



- {l2/x)V 
12 + 2x2 



Again, regularity at x = requires Di — 0. 

In the super-horizon, x <C 1 regime we obtain 



* = *o, D 



with 



D, 







g-Do- -Vox\ V = Vox 



-6*0 = -D2 



(3.20) 
(3.21) 
(3.22) 



(3.23) 

(3.24) 
(3.25) 



(3.26) 

(3.27) 
(3.28) 



On sub-horizon, x ^ 1 scales we find oscillating solutions with constant amplitude 
with a frequency of 1/ VS: 



V = y,sin( ^ 

X 

Dn = Do COS ( 



D2^ 



V3 



(3.29) 
(3.30) 
(3.31) 
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Note that also for radiation perturbations 

D = ~VoX^ < * 

is small on super horizon scales, a; -C 1. The perturbation amplitude is given 
by the largest gauge invariant perturbation variable. We conclude therefore that 
perturbations outside the Hubble horizon are frozen to first order. Once they enter 

the horizon they start to collapse, but pressure resists the gravitational force and the 
radiation fluid starts to oscillate. The perturbations of the gravitational potential 
oscillate and decay like 1/a^ inside the horizon. 

3.3 Adiabatic and isocurvature initial conditions 
for a matter radiation fluid 

In this section we want to investigate a system with a matter and a radiation com- 
ponent that are coupled only by gravity. The matter component acts therefore as 
dark matter, since it does not interact directly with the radiation. 

Since the matter and radiation perturbations behave in the same way on super- 
horizon scales, 

D^;'^ + Bx^, ^A' + oc V^""^ oc X, (3.32) 

we may require a constant relation between matter and radiation perturbations. As 
we have seen in the previous section, inside the horizon {x > 1) radiation perturba- 
tions start to oscillate while matter perturbations keep following a power law. On 
sub-horizon scales a constant ratio can thus no longer be maintained. There are two 
interesting possibilities: 



3.3.1 Adiabatic initial conditions 

Adiabaticity requires that matter and radiation perturbations are initially in perfect 
thermal equilibrium. This implies that their velocity fields agree (see below, section 
on Boltzmann eqn.!) 

SO that the energy flux in the two fluids is coupled initially. 

Let us investigate the radiation solution in the matter dominated era, when 
the corresponding scale is already sub-horizon. Since ^ is dominated by the matter 
contribution, we have ^ ~ const = ^'o- We neglect the (decaying) contribution from 
the sub-dominant radiation to Energy-momentum conservation for radiation 
then gives 

dW' = -lyW (3.34) 
3 

= 2* + ^L»W. (3.35) 
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Now \1/ is just a constant given by the matter perturbations, and it acts like a 
constant source term. The full solution of this system is then 



X \ 




^ X 






X \ 


+ — Asm 


( X 


7!; 





V^M = 5 cos 
Our adiabatic initial conditions require 



- 2\/3^' sin 



X 



(3.36) 
(3.37) 



y(r) y(m) 

lim = Vo = lim < oo. 



x^O X 



x-^O X 



(3.38) 



Therefore 5 = and A = 4Vo - 8^. Using in addition ^ = Wq (see (|]T|)) we 
obtain 



yir) 

9 



44^ 

w cos 

3 



(- 



X 

-^(5 + ^x2) 



3 

31^0- 



8^ 



On super-horizon scales, x -C 1 we have 



20 



D\'^ ^ m and V^'^ 

^ 3 



-x^ 



(3.39) 
(3.40) 

(3.41) 

(3.42) 
(3.43) 

(3.44) 



note that Z)^*^^ = (4/3) Dg™''' and V^^*"^ = K*^™"-' for adiabatic initial conditions. 



, (m) 



3.3.2 Isocurvature initial conditions 

Here we want to solve the system ( p.50| ) and ( |2.55|) for dark matter and radiation 
under the condition that the metric perturbations vanish initially, i.e., \1/ = 0, 



^ = — 



k- 



Dg + ?,{! + w)^! + ?,{! + w) 



0. 



(3.45) 



In principle, we have four evolution and one constraint equations. We therefore 
have four constants to adjust. Condition ( 3.45 ), however, requires an entire function 
to vanish. This may be impossible. Let us nevertheless try: 
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If ^ = the solutions of the radiation dominated equations are simply 

Z,.. . ^.„.(^)^B.„(^) (3.46) 

For the matter perturbations we find 

y{m) ^ _Y^^ a(xx^, 1<P<2 (3.48) 

a 

^ p — la 

j^{ra) ^ C^"') - Vo\og{x) if/3 = l (3.50) 

Here (3 is the exponent of the scale factor a (x r]^, hence /3 = 1 in the radiation era 
and /? = 2 in the matter era. 
^ = implies with 

Dg = ^(p,L>f +p^L>f)) and (3.51) 

V = -^{{Pr+Pr)V^'-^+PmV^"'^) (3.52) 
P+P 



that 



(3.53) 



Since V^""^ ocl/ait can compensate, for small values of x, the term oc cos(x/ y/S) 
of V^'^\ which behaves like 1/a as well, due to the pre-factor Pr/Pm- This term can 
also be compensated in Dg ^ by the term Vox/a of Dg™"^ In the purely radiation dom- 
inated universe, the log-dependence of Dg"^^ renders this compensation imperfect. 
However, there is no way to compensate C'-"*-' or the term proportional to A. We 
therefore have to choose A = C^"*^ = and 

afL}^B^Vo, B=P^VW,. (3.54) 

Pm 3 apr 

(The compensation of the smaller terms in D^g'' and Dg"*"* is only complete if /3 ~ 2.) 
With Cs = 1/^3 we find 

sin [cgx) (isocurvature) (3.55) 

aprCs 

~ * ^8 - y COS {csx)^ (adiabatic) . (3.56) 
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The CMB anisotropies, which we are going to determine in the next chapter, contain 
a term 

AT 1 

— (k, no,n) = --- + (k, r^dec) e^"^-^^'-^--^ ■■■ . (3.57) 

On scales where this term dominates, the peaks in Dg translate into peaks in the 
angular power spectrum of CMB anisotropies. 

For isocurvature initial conditions, we find a first peak in Dg at 

Here 'df'^ is the angle under which the comoving scale A^^"* at comoving distance 
Vo "~ ^dec is seen. In the next chapter, we will expand the temperature fluctuations 
in terms of spherical harmonics. An fluctuation on the angular scale then shows 
up around the harmonic £ ~ vr/^?. As an indication, we note that for A — k — 0, 
the harmonic of the first isocurvature peak is 

41) ^ ~ 110 , 

In the adiabatic case the first "peak" is at k'^a^ = 0. 

Since D^g"^ is negative for small x, the first peaks are "expansion peaks", and 
due to the gravitational attraction of the baryons (which we have neglected in this 
simple argument) they are less pronounced than the second ( "compression" ) peaks. 

These second peaks are usually called the "first acoustic peak". (It is the first 
compression peak and we shall adopt the convention to call it the "first peak" mainly 
for consistency with the literature.) They correspond to wavelengths and angular 
scales 

Af ^ = -c,77dec , 4^ - ^y^^^-^^- , ^ 350 (isocurvature) (3.59) 

3 X (^0 - Vdec) 

Ai') = c,77dec, ^(2) ^-^'^-^ , £(2)^ 220 (adiabatic). (3.60) 

X [VO - Vdec) 

Here the indicated harmonic is the one obtained in the case A = k = 0, for a typical 
baryon density inferred from nucleosynthesis. 

It is interesting to note that the distance between consecutive peaks is indepen- 
dent of the initial condition. It is given by 

Aki^kl'^-kl'^=7i/{csVdec)=Aka, At? = — ^f!?^^ , A£ ~ 220 . (3.61) 

X [m - Vdec) 

Again, the numerical value indicated for Ai corresponds to a universe with A = 
K — O. The result is strongly dependent especially on k. This is the reason why the 
measurement of the peak position (or better of the inter-peak distance) allows an 
accurate determination of curvature. 
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From our analysis we can draw the following important conclusions: For scales 
where the Dg^''-term dominates, the CMB anisotropies show a series of acoustic 

(2) 

oscillations with spacing Ak, the position of the first significant peaks is at A; = k^^^, 
depending on the initial condition. 

The spacing Ak is independent of initial conditions. The angle Ai} onto which 
this scale is projected in the sky is determined entirely by the matter content and 
the geometry of the universe. According to our findings in Chapter I, will be 
larger ii fl^ < (positive curvature) and smaller if fl^ > (see Fig. |L3l) . 

In our analysis we have neglected the presence of baryons, in order to obtain sim- 
ple analytical results. Baryons have two effects: They lead to (p+3p)rad+bar > 0, and 
therefore to an enhancement of the compression peaks (the first, third, etc. acous- 
tic peak). In addition, the baryons slightly decrease the sound speed c^, increasing 
thereby Ak and decreasing A-d. 

Another point which we have neglected is the fact that the universe became 
matter dominated at rjeq, only shortly before decouphng: ?7dcc — 4r7cq, for flm = 1- 
As we have seen, the gravitational potential on sub-horizon scales is decaying in 
the radiation dominated era. If the radiation dominated era is not very long before 
decoupling, the gravitational potential is still slightly decaying and free streaming 
photons fall into a deeper gravitational potential than they have to climb out of. This 
effect, called "early integrated Sachs- Wolfe effect" adds to the photon temperature 
fluctuations at scales which are only slightly larger than the position of the first 
acoustic peak for adiabatic perturbations. It therefore 'boosts' this peak and, at the 
same time, moves it to lightly larger scales (smaller angles). Since r^gq oc /i"^, the 
first acoustic peak is larger if h is smaller. 

A small Hubble parameter increases therefore the acoustic peaks. A similar effect 
is observed if a cosmological constant or negative curvature are present, since r/eq is 
retarded in those cases. 

The real universe contains not only photons and dark matter, but also neutrinos 
and baryons. It has actually be found recently |T6| that this 4 fluid mixture allows 
five different modes which grow or stay constant. The adiabatic mode, the dark 
matter isocurvature mode which we have just discussed, a similar baryon isocurva- 
ture mode and two neutrino isocurvature modes. The most generic initial conditions 
which allow for arbitrary correlations between the different modes are very unpre- 
dictable. We can maybe just say that they lead to a first acoustic peak in the 
range of 150 < i^'^^ < 350 for a spatially fiat universe. In the rest of this review, 
we only discuss adiabatic perturbations, which are by far the most studied, but 
it is important to keep in mind that all the results especially concerning the esti- 
mation of cosmological parameters is not valid if we allow for more generic initial 
conditions |^ . 
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3.3.3 Vector perturbations of perfect fluids 

If n(^) = equation (|]5|) implies 

Q oc a^''-\ (3.62) 

For p/p = < Vs) this leads to a non-growing vorticity. The dynamical Einstein 
equation implies 

cr(^) oc a"^ , (3.63) 
and the constraint ( |2.51| ) reads (at early times, so we can neglect curvature) 



(3.64) 

If perturbations are created in the very early universe on super-horizon scales 
{e.g. during an inflationary period), vector perturbations of the metric decay and 
become soon entirely negligible. Even if f2j remains constant in a radiation domi- 
nated universe, it has to be so small on relevant scales at formation {xin ^ 1) that 
we may safely neglect it. 

3.3.4 Tensor perturbations 

The situation is different for tensor perturbations. Again we consider the perfect 



fluid case, II-J^ = 0. There ( |2.54| ) implies (if k is negligible) 



^ + ^H[^ + Hij=0, (3.65) 

with /3 = 1 in the radiation dominated era and /5 = 2 in the matter dominated era. 
The less decaying mode solution to Eq. (|3.65|) is Hij = eijX^^'^~'^Ji/2-p{x), where 
denotes the Bessel function of order z/ and e^j is a transverse traceless polarization 
tensor. This leads to 

Hij = const for X ^ 1 (3.66) 

Hii = - for X ~ 1. (3.67) 
a 
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Chapter 4 

CMB anisotropies 



4.1 Light-like geodesies 

After decoupling, rj > rjdcc, photons follow to a good approximation light-like geodesies. 
The temperature shift is then given by the energy shift of a given photon. 

The unperturbed photon trajectory follows (x^) = {rj, n(?7 — 770) +Xo), where xq is 
the photon position at time r/o and n is the (parallel transported) photon direction. 
With respect to a geodesic basis (e)^^^, the components of n are constant. U k — 
we may choose ej = d/dx'^; if k 7^ these vector fields are no longer parallel 
transported and therefore do not form a geodesic basis (Ve^ej = 0). 

Our metric is of the form 

ds^ = a'^ds'^ ,with (4.1) 
ds'^ = (7^^ + hf,^) dx^dx", 7oo = -1, 7^0 = 0, 7^^ = 7^^ (4.2) 

as before. 

We make use of the fact that light-like geodesies are conformally invariant. More 
precisely ds"^ and ds"^ have the same light-like geodesies, only the corresponding 
affine parameters are different. Let us denote the two affine parameters by A and A 
respectively, and the tangent vectors to the geodesic by 

n = ^, n = S , = = , n° = 1 , = 1. (4.3) 
dX dX 

We set 71^ — 1 + SnP. The geodesic equation for the perturbed metric 

ds'' = (7^, + h^,)dx^dx'' (4.4) 

yields, to first order, 

A^^M = _(5r^ (4.5) 
dX 

For the energy shift, we have to determine 5n^. Since g^^ = — 1 • 5o/i + first order, 
we obtain 5r°^ = - V2(/iao|/3 + /i/soia - ^a/?), so that 

j^5n' = K,\pn^n'' - ^Kpn^n^ . (4.6) 
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Integrating this equation we use hao\i3n^ = -^{haon"), so that the change of n° 
between some initial time rji and some final time i]f is given by 



Sn^\{ = [hoo + ho,n^]{ -\ f h.^n^n^dX . (4.7) 



On the other hand, the ratio of the energy of a photon measured by some observer 
at t/ to the energy emitted at ti is 

Ef _ {n- u)f _Tf {n- u)f ^^^^ 
Ei (n • u)i Ti (n ■ u)i ' 

where Uf and Ui are the four- velocities of the observer and emitter respectively, and 
the factor T//Tj is the usual (unperturbed) redshift, which relates n and n. The 
velocity field of observer and emitter is given by 

u = (1 - A)dr, + v'di . (4.9) 

An observer measuring a temperature Tq receives photons that were emitted 
at the time rjdec of decoupling of matter and radiation, at the fixed temperature 
Tdec- In first-order perturbation theory, we find the following relation between the 
unperturbed temperatures Tf, Ti, the measurable temperatures Tq, Tdec, and the 
photon density perturbation: 

where 6^^^ is the intrinsic density perturbation in the radiation and we used 
in the last equality. Inserting the above equation and Eq. ( [4. 71 ) into Eq. (|4.8|), and 
using Eq. ( p.l9| ) for the definition of h^^^, one finds, after integration by parts the 
following result for scalar perturbations: 



^ To_ 
Ei Td^Q 



-D^;^ + Vj^^n^ + ^ - $ 



+ [\^-^)dX} . (4.11) 



Here Dg^^ denotes the density perturbation in the radiation fiuid, and V^'^'' is the 
peculiar velocity of the baryonic matter component (the emitter and observer of 
radiation). The final time values in the square bracket of Eq. ( [4.11|) give rise only 



to monopole contributions and to the dipole due to our motion with respect to the 
CMB, and will be neglected in what follows. 

Evaluating Eq. ( |4.11| ) at final time tjq (today) and initial time ridec, we obtain 



the temperature difference of photons coming from different directions n and n' 

^ ^ «^ _ ,4.12) 
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with temperature perturbation 



AT(n) 
T 



D^:^ + V}^^n^ + ^ - $ 



rvo 

iVdec, Xrfec) + / - ^){v,Av))dv , (4.13) 

•^Vdec 



where x(?7) is the unperturbed photon position at time t] for an observer at Xq, 
and Xdec = ^iVdec) (If K = we simply have x(?7) = xq — (?7o — ''7)n.). h The 
first term in Eq. ( [4. 131 ) describes the intrinsic inhomogeneities on the surface of 
last scattering, due to acoustic oscillations prior to decoupling. Depending on the 
initial conditions, it can contribute significantly on super-horizon scales. This is 
especially important in the case of adiabatic initial conditions. As we have seen in 
Eq. ( |3.44| ), in a dust + radiation universe with f2 = 1, adiabatic initial conditions 
imply D^g\k,r]) = -20/3^(A;, r/) and 1/^^) = V'-''^ < for kr] < 1. With $ = 
the the square bracket of Eq. ([4.13|) gives 



-W(?7rfec,Xrfec) 



T 3 

■' adiabatic 



on super-horizon scales. The contribution to ^ from the last scattering surface 
on very large scales is called the 'ordinary Sachs- Wolfe effect' (OSW). It has been 



derived for the first time by Sachs and Wolfe For isocurvature perturbations, 
the initial condition D^g \k, rj) ^ for rj ^ is satisfied and the contribution of Dg 
to the ordinary Sachs- Wolfe effect can be neglected. 



\ / isocurvature 



The second term in (|4.13|) describes the relative motions of emitter and observer. 
This is the Doppler contribution to the CMB anisotropics. It appears on the same 
angular scales as the acoustic term, and we thus call the sum of the acoustic and 
Doppler contributions "acoustic peaks". 

The last two terms are due to the inhomogeneities in the spacetime geometry; the 
first contribution determines the change in the photon energy due to the difference 
of the gravitational potential at the position of emitter and observer. Together 
with the part contained in dI['' they represent the "ordinary" Sachs- Wolfe effect. 
The integral accounts for red-shift or blue-shift caused by the time dependence of 
the gravitational field along the path of the photon, and represents the so-called 
integrated Sachs- Wolfe (ISW) effect. In a = 1, pure dust universe, the Bardeen 
potentials are constant and there is no integrated Sachs- Wolfe effect; the blue-shift 
which the photons acquire by falling into a gravitational potential is exactly canceled 
by the redshift induced by climbing out of it. This is no longer true in a universe 
with substantial radiation contribution, curvature or a cosmological constant. 

The sum of the ordinary Sachs-Wolfe term and the integral is the full Sachs- 
Wolfe contribution (SW). 
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For vector perturbations S^^^ and A vanish and Eq. ( [4.8|) leads to 



(4.14) 



We obtain a Doppler term and a gravitational contribution. For tensor perturba- 
tions, i.e. gravitational waves, only the gravitational part remains: 



rf 

{Ef/E,Y^^ = (a,/a/)[l - / Hi.n'n^dX] . 

J i 



(4.15) 



Equations ( ^.11| ), ( [4.14[ ) and ( [4.15[ ) are the manifestly gauge invariant results for the 
Sachs-Wolfe effect for scalar vector and tensor perturbations. Disregarding again 
the dipole contribution due to our proper motion, Eqs. ( [4.14| , |4.15| ) imply the vector 
and tensor temperature fluctuations 



ATfnl 



(V) 



V/"'' (?7dec, ^dec)n^ + / CTj {ill ^{ri))n^d\ 



T 

/ AT(n) ^ 



(4.16) 
(4.17) 



Note that for models where initial fluctuations have been led down in the very early 
universe, vector perturbations are irrelevant as we have already pointed out. In this 
sense Eq. ( [4.16| ) is here mainly for completeness. However, in models where pertur- 
bations are sourced by some inherently inhomogeneous component [e.g. topological 
defects) vector perturbation can be important. 



4.2 Power spectra 

One of the basic tools to compare models of large scale structure with observations 
are power spectra. They are the "harmonic transforms" of the two point correla- 
tion functions. If the perturbations of the model under consideration are Gaussian 
(a relatively generic prediction from inflationary models), then the power spectra 
contain the full statistical information of the model. 

One important power spectrum is the dark matter power spectrum, 

P,,(fc) = (|Df)(k,r/o)f), (4.18) 

where ( ) indicates a statistical average over "initial conditions" in a given model. 
Pd(/c) is usually compared with the observed power spectrum of the galaxy distri- 
bution. 

Another power spectrum is given by the velocity perturbations, 

Pv{k) = (|V(k,r7o)|'> ^ H',Q'-'PDik)k~' . (4.19) 
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For ~ we have used that |fc\^|(?7o) = D^g^\r]Q) ~ HoVt^'^Dg on sub-horizon scales 
(see e.g. Q). 

The power spectrum we are most interested in is the CMB anisotropy power 
spectrum. It is defined as follows: AT/T is a function of position xq, time rjo and 
photon direction n. We develop the n-dependence in terms of spherical harmonics. 
We will suppress the argument tiq and often also xq in the following calculations. 
All results are for today (rjo) and here (xq). By statistical homogeneity expectation 
values are supposed to be independent of position. Furthermore, we assume that 
the process generating the initial perturbations is statistically isotropic. Then, the 
off-diagonal correlators of the expansion coefficients vanish and we have 

AT 

— (xo, n, rjo) = ^ a£„(xo)Y£m(n), {a^m ■ a\im) = ^u'^mm'Ci (4.20) 

£,m 

The Cis are the CMB power spectrum. We assume that the perturbations are 
generated by a homogeneous and isotropic process, so that Ci depends neither on 
Xo nor on m, and that {a£m ■ ci£'m') vanishes for £ 7^ or m 7^ m'. 

Let us, at this point insert a comment on the problem of cosmic variance: Even 
if our 'ergodic hypothesis' is correct and we may interchange ensemble and spatial 
averages, we cannot obtain very precise averages for measurements of large scale 
characteristics, due to the fact that we can observe only the universe around a given 
position. For example, let us assume that temperature fiuctuations are Gaussian, 
as they are in most inflationary models. The functions aim are then also Gaussian 
distributed, and we have a variance of 



1 ^ 
2i+l ^ 



|2 



-a 



c 

10^ — 0/> 



2£+ 1 



on the average of the 2£ + 1 values which can in principle be measured from one 
point with full sky coverage. For simplicity, we neglect the increase of the variance 
due to the fact that our own milky ways blocks a portion of sky of about 20%. 
Wick's theorem now gives 

For a given multipole I we then expect a variance of the C^'s 

Ci = V^' ^^-^^^ 

in real experiments, this 'cosmic variance' is in general much larger due to the limited 
sky coverage. 

The two point correlation function is related to the C^'s by 
/ AT AT \ 

^— (n) — (n')^ = ^ (a£„-a^,„,)y<.^(n)F;„(n') = 
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^ 1 

Y^QY, y,„(n)y4(n') = Yi^i+l)CePe{fi), (4.22) 



An 

m=-£ e 



2|±iP.(n-n') 



where we have used the addition theorem of spherical harmonics for the last equality. 
The PiS are the Legendre polynomials. 

Clearly the a^m's from scalar, vector and tensor perturbations are uncorrelated, 

(4i'4;2.) = = ("Z!"?:,.) = o. (4.23) 

Since vector perturbations decay, their contributions, the C^^\ are negligible in 
models where initial perturbations have been laid down very early, e.g., after an 
inflationary period. Tensor perturbations are constant on super-horizon scales and 
perform damped oscillations once they enter the horizon. 

Let us first discuss in somewhat more detail scalar perturbations. We specialize 
to the case k = for simplicity. We suppose the initial perturbations to be given by 
a spectrum, 

= A^fc"-^//^-^ (4.24) 

We multiply by the constant ?7o~^, the actual comoving size of the horizon, in order 
to keep A dimensionless for all values of n. A then represents the amplitude of 
metric perturbations at horizon scale today, k = l/rjQ. 

On super-horizon scales we have, for adiabatic perturbations: 

i^M = -^^ + 0{x^), V^'^ = V^'^ = 0{x) (4.25) 
4 3 

The dominant contribution on super-horizon scales (neglecting the integrated 
Sachs- Wolfe effect J $ - ^ ) is then 

AT 1 

— (xo, n, r]o) = -^'(xdec, Vdec)- (4.26) 
The Fourier transform of ( [4.26| ) gives 



AT 



Using the decomposition 



(k, n, r]o) = l^ik, r^dec) ■ e^'^^^""""'^-) . (4.27) 



3ikn(j?o-f?dcc) _ 



^(2£ + l)t'je{k{r]o - Vdcc))Pi{k ■ n) 
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where je are the spherical Bessel functions, we obtain 
/AT AT \ 

( — (xo, n, r^o) — (xo, n , r]o) ) (4.28) 

1 f . I AT AT \ 

= y c^'^xo ( — (xo,n,r/o) — (xo,n',r/o)^ 

IP °° 
^ e,e'=o 

■Pi{k ■ n) ■ P;(k • n') . (4.29) 

Inserting P,(kn) = ^.EmY^mimM andP;(kn') = ^ y;^,(k)y,,^,(n'), 
integration over the directions dQf^ gives 5tt'5mm' Xlm ^^m(^)^m('^')- Using as well 
^^y/^(n)F^„(n') = ^P^(yLi), where /i = n ■ n', we find 



AT AT 

— (xo, n, r/o) — (xo, n', r]o) 



nn'=/i 



v:^2^+l , ,2 /"c/A; /I, ,,\ , 



Comparing this equation with Eq. (|4.22|) we obtain for adiabatic perturbations 



on scales 2 < £ < xivo - Vdcc)hdec ~ 100 



3 



2\ 



k^'jUk iVo-Vdcc))- (4.31) 



If \1' is a pure power law and we set k{r]Q — rj^ec) ~ krjQ, the integral ( [4.31 
be performed analytically. For the ansatz ( [4.24|) one finds for —3 < n < 3 



can 



risw) ^ A' r(3-n)r(£-i + t) 
^ 9 23— r2(2-f)r(£ + f-f)' ^ ■ ' 

Of special interest is the scale invariant spectrum, n = 1. This spectrum with 
a time and scale independent gravitational potential has first been investigated by 
Harrison and by Zel'dovich p5|. It is therefore called the Harrison-Zel'dovich spec- 
trum. It leads to 

i{i + l)Cf ^) = const. :^ (^0^^)) ^ ' ^i^n/£. (4.33) 

This is precisely (within the accuracy of the experiment) the behavior observed by 
the DMR experiment aboard COBE g]. 
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Inflationary models predict very generically a HZ spectrum (up to small cor- 
rections). The DMR discovery has therefore been regarded as a great success, if 
not a proof, of inflation. There are however other models like topological defects 



p8| , p9|] or certain string cosmology models I^O] which also predict scale-invariant, 
i.e., Harrison Zel'dovich spectra of fluctuations. These models do however not be- 
long to the class investigated here, since in these models perturbations are induced 
by seeds which evolve non-linear ly in time. 

For isocurvature perturbations, the main contribution on large scales comes from 
the integrated Sachs-Wolfe effect and ( [4.31| ) is replaced by 



(ISW) 



8 fdk_ 



'^ik,ri)j]{k{rio~ri))dri 



Vdc 



(4.34) 



Inside the horizon is roughly constant (matter dominated). Using the ansatz 
( [4.24| ) for \E' inside the horizon and setting the integral in (|4.34| ) ~ 2'${k,T] = 
l/k)j]{kriQ), we obtain again (|4.32| ), but with A^/d replaced by AA^. The Sachs- 
Wolfe temperature anisotropics coming from isocurvature perturbations are there- 
fore about a factor of 6 times larger than those coming from adiabatic perturbations. 

On smaller scales, i ~ 100 the contribution to AT/T is usually dominated by 
acoustic oscillations, the first two terms in Eq. (|4.13| ). Instead of ( [4.34|) we then 
obtain 



2 

./O 



dk 
T 



\D^\Kr]^,:)3t{kr]^) + V^'\Kr]^,,)3\{kr]^] 



4 



(4.35) 



On sub-horizon scales D^p and V^*^^-* are oscillating like sine or cosine waves 
depending on the initial conditions. Correspondingly the C^^*^^ will show peaks and 
minima. On very small scales they are damped by the photon diffusion which takes 
place during the recombination process (see next section). 

For gravitational waves (tensor fluctuations), a formula analogous to ( |4.32| ) can 
be derived (see appendix). 



IT 



Vd 



dvH{v,k)— — 

{k{r]o-v)r 



(f + 2)! 
(£-2)!' 



(4.36) 



To a very crude approximation we may assume if = on super-horizon scales 
and / driHji{k{r]o — t])) ~ ii(?7 = l/k)ji{kr]o). For a pure power law, 

|2\ 



this gives 



(T) 



kmH{k,7^=i/k)f)=Aik^-7^~-\ 

^ 2 (1 + 2)1 [dx 3j{x) 
- 7r(£-2)!^./ X -4 



(4.37) 



X'* 
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Figure 4.1: A COBE normalized sample adiabatic( solid line) and isocurvature 
(dashed line) CMB anisotropy spectrum, C.{C. + 1)C£, are shown on the top panel. 
The quantity shown in the bottom panel is the ratio of temperature fluctuations for 
fixed value of A (from Kanazawa et al. pl|). 

^ (^ + 2)! ,2 r(6-n^)r(f-2 + ^) 

(£-2)! ^26-"Tr2(Z-nr)r(£ + 4-^)' ^' ' 

For a scale invariant spectrum (n-p = 0) this results in 

The singularity at £ = 2 in this crude approximation is not real, but there is some 
enhancement of £(£ + l)Cf'^ at £ ~ 2. 

Since tensor perturbations decay on sub-horizon scales, £ ~ 60, they are not very 
sensitive to cosmological parameters. 

Again, inflationary models (and topological defects) predict a scale invariant 
spectrum of tensor fluctuations (jit ~ 0). 

On very small angular scales, ^ ~ 800, fluctuations are damped by coUisional 
damping (Silk damping). This effect has to be discussed with the Boltzmann equa- 
tion for photons derived in the next section. 



4.3 The Boltzmann equation 
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Scalar 



Tensor 




Figure 4.2: Adiabatic scalar and tensor CMB anisotropy spectra are shown (top 
panels). The bottom panels show the corresponding polarization spectra (see Sec- 
tion I 



(from pi). 



4.3.1 Elements of the derivation 

When particles are not very tightly coupled, the fluid approximation breaks down 
and they have to be described by a Boltzmann equation, 



P'd.f - KpP'-p'^ = c[f] 



(4.40) 



C[f] is a collision integral which describes the interactions with other matter compo- 
nents. The left hand side of (|4.40|) just requires the particles to move along geodesies 
in the absence of collisions. 

Let us flrst consider the situation where collisions are negligible, C[f] = 0. The 
unperturbed Boltzmann equation implies that / be a function of v = ap only. Setting 
/ = f{v) + F(?7, X, f, n), where n denotes the momentum directions, leads then to 
the perturbation equation 



Bi\j Hij 



Hr 



(4.41) 



Here F^'^''* are the Christoffel symbols of the space of constant curvature k. 



To derive ( [4.41|) , we have used = 0. The Liouville equation for particles with 
non- vanishing mass can be found in Ref. 0. 
The ansatz 

Tv 



fix,p) = f 



^(^)(p,p)2 

T(x, n) 



/ 



T{x, n) 



(4.42) 



38 



with T(x, n) = T{ri) + AT(x, n) leads to 



f = f-v 



dfAT 



dv T 

Integrating ( [4.42| ) over photon energies, we can also write 

AT _ 1 

4'' 

where t is the brightness perturbation, 

477 

I = —r I Fv'^dv. 



pa^ Jo 



Comparing Eq. ([4.43|) with (|4.41| ), we find 



\ T 



(4.43) 



(4.44) 



(4.45) 



(4.46) 

The fact that gravitational perturbations of Liouville's equation can be cast 
purely in temperature perturbations of the original distribution is not astonishing. 
This is just an expression of gravity being "achromatic", i.e. independent of the 
photon energy. 

We now decompose ( ^.46| ) into scalar, vector and tensor components. Even 
though AT/T is just a function, it can be represented in the form 



AT 



X, n) = y^aii,...,i,(x)n'' 



■ n 



(4.47) 



where the are symmetric traceless tensor fields that contain scalar, vector, 

2-tensor and in principle also higher tensor components. Since spin components 
larger than 2 are not sourced by the right hand side of equation ( |4.46 ) and since 
they are suppressed at early times, when collisions are important, we neglect them 
here. 

For the scalar contribution to AT/T we obtain from ([4.46|) 



/AT\ 



(S) 



/AT\ 



(S) 



^{S)i j k 



T_ 

kfiA + fi'^k'^ [B - Ht] + Hl + -k^Hr 



where we have introduced the "direction cosine" /x defined by n^Y,i 
that in fiat space, k = 0, we have just fi = ik ■ n. 



(4.48) 
kfxY. Note 
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This equation is not manifestly gauge-invariant. However, setting 



A<(^)=f— J ^Hi^^-k^ET^ky.(H-BY (4.49) 

it reduces to 



dr^M^^^ + A;/xA^(^) - rg) Vn'^^^:^^ = ($ - ^) , (4.50) 

where $ and \1/ are the Bardeen potentials. If are components with respect to a 
geodesic basis (or k = 0), the third term on the left hand side of Eq. ( |4.50| ) vanishes. 
For simplicity we now concentrate on the case k = 0. We can then integrate the 
equation and obtain 

A^(^)(r7o,n,k) = exp[ik-n(r/i„-r/o)]A<^^H^m,n,k) 

rm 

+ / iexp[zk ■ n(?7 - r7o)]n ■ k ($ - \&) (i?7 . (4.51) 

Integration by parts and neglecting the monopole term ($ — \E') (r^o), leads to 

A<(^)(r7o,n,k) = 

exp[ik • Yi{r]in - r^o)] [M^^\riin, n, k) + ($ - ^) (?],„, k)] 

- / exp[?k ■ n(?7 - rjo)] - dr] . (4.52) 

Comparing this equation with ( [4.13|) , we see again that A^*^"^) = (■^)'"'^'* (up to 
gauge dependent monopole and dipole contributions) if the initial condition is 

M^'\v^n, n, k) = \D^;\v^n, k) + n ■ kV^'^V^r,, k) , 

which is equivalent to require that Ai^^\r]in) has no higher than first moments. As 
we shall see below, this assumption is quite reasonable since collisions suppress the 
build up of higher moments before recombination. 

Since the right hand side of ( [4.50| ) is gauge invariant, the left hand side must be 
so as well and we conclude that Ai^^^ is a gauge-invariant variable (a direct proof 
of this, analyzing the gauge transformation properties of the distribution function, 
can be found in Ref. 0). 

Ai^^^ used in this work coincides with the scalar temperature fluctuations up 
a to a gauge dependent monopole and dipole contribution. In other work the 
gauge invariant variable Q = Ai^^^ + ^ has been used. Since $ is independent of the 
photon direction n this difference in the definition shows up only in the monopole, 
Co. 
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The vector and tensor parts of 1ST jT are gauge-invariant by themselves and we 
denote them by A^*^^-* and M.^"^^ for consistency. In the absence of coUisions and 
with vanishing spatial curvature, they satisfy the equations 



• I m TT 

-m n Hmi- 



(4.53) 
(4.54) 



The components of the energy momentum tensor are obtained by integrating the 
second moments of the distribution function over the mass shell, 



p^dpdflp 



(4.55) 



'Pm{x) 

where Qp denotes the angular integration over momentum directions. One finds for 
K — and setting = n • k: 
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niUjM^^^dQ,. 



(4.56) 
(4.57) 

(4.58) 

(4.59) 
(4.60) 

(4.61) 

(4.62) 



Let us now turn to the collision term. At recombination (when the fluid de- 
scription of radiation breaks down) the temperature is ~ 0.4 eV. The electrons 
and nuclei are non-relativistic and the dominant collision process is non-relativistic 
Thomson scattering. Since collisions are important only before and during recombi- 
nation, where curvature effects are entirely negligible, we set k = in the reminder 
of this section. 

The collision term is given by 



C[f] 



dr) 



dr) ' 



(4.63) 



where /+ and /_ denote the distribution of photons scattered into respectively out 
of the beam due to Compton scattering. 
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In the matter (baryon/electron) rest frame, which we indicate by a prime, we 
know 

where Ue denotes the number density of free electrons, ctt is the Thomson cross 
section, and cu is the normahzed angular dependence of the Thomson cross section: 

3 1 
uj{n, n') = 3/4[l + (n • n'Y] = 1 + -ny-n-^- with riij = nirij — . 

In the baryon rest frame which moves with four velocity the photon energy is 
given by 

We denote by p the photon energy with respect to a tetrad adapted to the slicing 
of spacetime into {r] = constant} hyper-surfaces: 

p = ' with n = a-^[{l - A)dr, + B'di] . 

The unit vector n is the normal the the hyper-surfaces of constant time. The lapse 
function and the shift vector of the slicing are given hy a — a{l + A) and (3 — —B'^di 
. In first order, 

Po ^ ap{l + A) - apriiB^ , 

and in zeroth order, clearly. 

Pi = a,pni . 

Furthermore, the baryon four velocity has the form — a~^(l — A) , — nPv^ 
up to first order. This yields 

p'^p^u^^pii + miv'-B')). 

Using this identity and performing the integration over photon energies, we find 

dt' 4:71 J 

The distribution of photons scattered out of the beam, has the well known form 
(see e.g. Lifshitz and Pitajewski [1983]) 

df- f,( I X 

— = UTUef (p , n) , 

so that we finally obtain 

where 5r — (l/47r) J L{n)dD, is the photon energy density perturbation. 

Using the definitions of the gauge-invariant variables M.'^^^ and V'^^'^ for the photon 
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brightness perturbation and the baryon velocity potential, we can write C in gauge- 
invariant form. 



.1 



1 



(4.64) 



with dP = (I/tt) jM^^Un and 



8n I ^ ^ «J 



Since the term in square brackets of ( 4.64] ) is already first order we can set dt' = 
dt which yields C = = = aC. The Boltzmann equation for scalar 

•' drj dri ^ 

perturbations expressed in terms of the gauge invariant variable Ai^^'' then becomes 



M^'^^+n'diM^'^^ = n'di{<!>-<i/)+aaTne[]D^;^ -M^'^^ -n'diV^^^ + -n^jM'^] . (4.65) 

4 2 



For vector and tensor perturbations we obtain in the same way 



iv) 



n 



4.3.2 The tight coupUng Umit 

Before recombination, when He — Pb/rrip, 

1 10 



(4.66) 
(4.67) 



Vt 
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1 + Zeq) 2(1 + z) Z, 



eq • 



(4.68) 

(4.69) 



To lowest order in r/^, collisions force the photon distribution to be of the form 

(4.70) 



4 



the building up of higher moments is strongly suppressed by collisions. 

During recombination, the number density of free electrons, ne, decreases rapidly 
and the collision term becomes less and less important. Higher moments in the 
photon distribution develop by free streaming. 

The collision term C[A^'-'^''] of equation ( |4.65| ) also appears in the equation of 
motion of the baryons as a drag. The Thomson drag force is given by 



(4.71) 
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with Mj = — UjM^^'^dn . 
This yields the following scalar baryon equation of motion in an ionized plasma 

+ {d/a)V^'^ = k^- ^^^I^i-k^M, + V^'^) , (4.72) 

where we have added the drag force to the second eq. of ( p.55| ) with w = = 0. 

We now want to discuss equations ([4.65| , |4?7^ ) in the limit of very many collisions. 
The comoving photon mean free path is given by rj^ = It = {cLC^T^e)'^ ■ In lowest 
order ^t/v Zr/A, [| Al*-'^'' is given by (|4.7CI| ), and Eq. (|4.72| ) implies 

-^Mj + V^^^ = . 

Inserting the solution ( [4.70 ) in the Boltzmann equation ( [4.65 ) and integrating over 
directions this implies 

kV^b) ^ ^ j^y(r) ^ I^eXj) , (4.73) 

Implying especially V^^^ = V^^^ = V. Eq. ( [4.73[ ) is equivalent to the energy con- 
servation equation ([2.55[ ) for radiation. Using also ( |2.55[) for baryons, w = 0, we 
obtain 

^ 3 3 ^ ■ 

This shows that entropy per baryon is conserved, F = 0. Before recombination, 
when the collisions are sufficiently frequent, baryons and photons are adiabatically 
coupled. Inserting ([4.701) in ( [4.65[ ) we find up to first order in 7]t 

+in^kjD^;^ + An'n^kikjV + -n'rimjkiM"'^ - iAn^kj{^ - ^)] .(4.74) 

2 

Using ( [4.74[ ) to calculate the drag force yields 

F,=tk,ipr/3)[AV-D^;^+Ai^-n ■ 
Inserting Fj in ( [4.72[ ), we obtain 

{p, + {4/3)pr)V + pi>{d/a)V = {pr/S)D^;^ + {p, + (4/3)p,)^ - (4p,/3)<l' . 

This is equivalent to momentum conservation, the second equation of ([2.55[ ) for 
p = Pb + Pr, p = Pr/3 and r = n = 0, if we use 

= (4/3)D(^) and D„ = P''^' + ^^^^ . 

Pb + Pr 



^Here A is a typical size of a perturbation. For a given Fourier mode fc, it is A ~ ir/k. 
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In this hmit therefore, baryons and photons behave hke a single fluid with density 
P = Pr + Pb and pressure p = Pr/S. 

From (|2.55| ) we can derive a second order equation for Dg. This equation can be 
simplified if expressed in terms of the variable D related by ( |3.1| ). We obtain 



D+cik^D+{l+3ci-Qw){d/a)D-3[w{d/a)-{d/ay{3{ci-w)-{l/2){l+w))]D = . 

For small wavelengths (sub-horizon), which are however sufficiently large for the 
fluid approximation to be valid, l/r]T ^ Cgk ^ l/?7, we may drop the term in 
square brackets. The ansatz D{t) = A{t)exp{—i J kcgdt) then eliminates the term 
of order c^k'^. For the terms of order Cgk/rj we obtain the equation 

2 A/ A + (1 + 3c^ - Qw){d/a) + c./c, = . (4.75) 

For the case = w =const. , this equation is solved by A oc {kviY~^ with z/ = 
2/{3w + 1), i.e., the short wave limit. In our situation we have 



w 



Pr 



3(Pr + Pb) 

Pr (4/3)Pr 



and 



3(Pr + Pb) 4:pr + 3pb 

Cs/cs = -3/2(d/a)- — . 

Apr + 3pb 

Using all this, one finds that 

^_ , Pb+{mPr V _fp + P^'^' 



Cs{Pr + PbYa^J ycsp'^a'^ 

solves ( |4.75| ) exactly, so that we finally obtain the approximate solution for the, 
tightly coupled matter radiation fiuid on sub-horizon scales 

D{t) oc ' exp(-^A; j Csdr]) . (4.76) 

Note that this short wave approximation is only valid in the limit t] ^ l/{csk), thus 
the limit to the matter dominated regime, Cs — > cannot be performed. In the limit 
to the radiation dominated regime, —>■ 1/3 and p oc we recover the acoustic 
waves with constant amplitude which we have already found in the last subsection. 
But also in this limit, we still need matter to ensure 7]t = l/{aaTne) -C rj. In the 
opposite case, rjT 3> r], radiation does not behave like an ideal fiuid but it becomes 
collisionless and has to be treated with Liouville's equation (( |4.65| ) without the 
collision term). 
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4.3.3 Damping by photon diffusion 

In this subsection we discuss the Boltzmann equation in the next order, (r/T^/r/)^. In 
this order we will obtain the damping of fluctuations in an ionized plasma due to the 
flniteness of the mean free path; the non-perfect coupling. We follow the treatment 
by Peebles [|17| (using our gauge-invariant approach instead of synchronous gauge). 
Again we consider Eqs. ( [4.65| ) and ( [4.72| ), but since we are mainly interested in 
collisions which take place on time scales ?7r <^ ?7, we neglect gravitational effects 
and the time dependence of the coefficients. We can then look for solutions of the 
form 

V oc A4^^'^ oc exp{i{k • x — uji])) . 

In ( [4.65| ) and ( [4.72| ) this yields (neglecting also the angular dependence of Compton 
scattering, i.e., the term n^jM*-') 

A.f5) 1 Di'^^-iknV 

= ; r 4.77 

4 1 - tr]Tiuj - k • n) ^ ^ 

and 

riTkioV = {4pr/3pb){ikV + M) , (4.78) 

with M = (3/47r) / nM^^MVt. Integrating (|]B3|) over angles, one obtains + 
(i/3)kjM^ = 0. With our ansatz therefore k • M = 3ujDg^\ Using this after scalar 
multiplication of (|4.78| ) with k, we find, setting R = 3pb/4:pr, 



V 



Inserting this result for V in ( [4.77| ) leads to 



4 1 — irixiuj — kp) 



where we have set /x = k ■ n. This is the result of Peebles ^7j, where the same 
calculation is performed in synchronous gauge. Like there (§92), one obtains in 
lowest non-vanishing order ujt]t the following dispersion relation: Using 

1 fi nW 1 /•! 1 + ^^^/^ 

' M^'^dp=^, which yields 1 = - / ^'^"^"f dp 



2 7„i ^ 4 ' 2 7_i 1 - iririuj - kp) 

one finds 

co = iuo-il with a;o = fc/[3(l + i?)]i/2 ^ = (fc2^^/6):^_L5i^^±l) . (4.79) 

In the baryon dominated regime, R> 1, therefore 

7 ^ k'^riT/6 . (4.80) 
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(If the angular dependence of Thompson scattering is not neglected, the term |(-R + 
1) in Eq. ( [4 .791 ) becomes |(-R + 1). If also polarization is taken into account, one 
obtains ^{R + l).) 

Posing kdampVT/G = 1, this leads to a damping scale Adamp ~ VTiVdec)/'^, which 
is projected in the microwave sky to an angle 

,9 . . vAVdcc) 

"damp 



2X(?70 - Vdec) ' 

For K = this corresponds to a few arc minutes and to the harmonic number 

4amp = vr/^damp ^ t^tt^tH ^ ^^\y^' ^>,h . (4.81) 
20?7T(77dcc) 20 

This estimate is very crude since we are using the the approximation for rjT from 
the tight coupling regime just where coupling stops to be tight, and we assume 
an arbitrary value of rif, ~ O.lrib at the moment of decoupling. Both these errors 
enhance the value of i'damp somewhat. Numerical results give 

4amp ~ 800 - 1000 

in a ft = universe. In open (closed) universes, this scale (which of course also 
depends on flf,) is moved to larger (lower) i due to projection. A reasonable approx- 
imation for the damping harmonic is 



"damp 



1000 



0.02(1 -fi«)i/2 



Temperature fluctuations on smaller scales, i > ^damp are exponentially damped by 
photon diffusion. 



4.4 Polarization and moment expansion 

Thomson scattering is not isotropic. And what is more, for a non-isotropic photon 
distribution it depends on the polarization: Even if the incident photon beam is 
unpolarized, the scattered beam will be, unless the incident distribution is perfectly 
isotropic. In the previous section we have neglected this effect by summing over 
initial polarizations and averaging over final polarizations. Now we want to derive 
the difference in the Boltzmann equation taking into account polarization. 
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Polarization is usually characterized by means of the Stokes parameters Jig, |T^, 
For a harmonic electromagnetic wave with electric field 

E(x, t) = [siEi + £2^2) e'"(""-*) , (4.82) 
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where n, Si and £2 form an orthonormal basis and the complex field amplitudes are 
parameterized as Ej = aje^^\ the Stokes parameters are given by 



/ 

Q 
u 

V 



2 2 



— a2 

2aia2 cos(52 — Si] 
2aia2 sin(52 — ^i) 



(4.83) 
(4.84) 
(4.85) 
(4.86) 



/ is the intensity of the wave (whose perturbation t has been introduced in the pre- 
vious section), while Q is a measure of the strength of linear polarization (the ratio 
of the principal axis of the polarization ellipse). V measures circular polarization 
which is not generated by Thomson scattering and therefore V vanishes if the ini- 
tial circular polarization vanishes (which we assume). U is then determined via the 
identity P = + U^. For scalar perturbations also U vanishes. 

Since Q vanishes in the background, to first order it obeys the unperturbed 
Boltzmann equation. 



{S)i 



n-'n 



(4.87) 



where C is the collision integral. The same type of equation, with a somewhat 
different collision integral is satisfied by U . The collision integral for V does not 
couple to I,Q or U and hence \^ = is a consistent solution. 

An explicit derivation of the following Boltzmann hierarchy including polariza- 
tion is presented in Appendix ^. Here we just repeat the necessary definitions and 
the results. 

The brightness anisotropy Ai and the non-vanishing Stokes parameters Q and 
U can be expanded as 

2 

M{v,k,n) = ^ ^ >lf)(r/,fc)oGr(n), (4.88) 

£ m=-2 
2 

Qiv, k, n) ± ^U{v, k, n) = H (^^^ ± ^i?f^)2Gr(n). (4.89) 



m=-2 



The special functions are described in Appendix 0. The coefficients m = 0,m = 
±1 and m = ±2 describe the scalar (S), vector (V^) and tensor (T) components 
respectively. The Boltzmann equation for the coefficients xj™'' is given by 



M'r'' - k 



2i 



2£ + 3 
{£ > m) 



(4.90) 



E 



(m) 



2i- 



-E 



2m 



(m) 



2^ 

21- 



-1 7-,(m) 
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B^r'^ - k 



2l^l p(m) 



2£- 1 



p(m) _ 2l^T+l p(m) 



(4.91) 



(4.92) 



where we set 



,(0) 



?(0) 



5, 



(0) 



n, 



5« = nearaC(« + 4S, 5. 



Sf^ = ne(JTCi4:UJb, 



(4.93) 



(2) 



n, 



.araC^^^ + 4ij 



and C^"*) = 7k[A^2"'^ - Ve^^^""^]. The couphng coefficients are 



(£2-m2)(£2-s2) 
£2 • 



The CMB temperature and polarization power spectra are given in terms of the 
expansion coefficients m!(^\ and 5^™-* as 



(2£ + l)2cf''(™) = / k'dkX^'^^Yl"'^* 

OTT J 



(4.94) 



where = 1 for m = and rim — 2 ior m — 1,2, accounting for the number of 
modes. Since B is parity odd, the only non-vanishing cross correlation spectrum is 



10- 

10- 

10- 
10- 
B 10- 

- 10- 

10- 
10- 

10- 
10- 




Figure 4.3: The temperature anisotropy (solid), the polarization (dashed) and their 
correlation (dotted) are shown for a purely scalar standard CMD model. 
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The Boltzmann hierarchy presented above can be solved numerically with pub- 
licly available fast numerical codes like CMBfast g or CAMCODE g. This 
enables us to compute the CMB anisotropy and polarization spectra for many dif- 
ferent values of cosmological parameters, and compare them with present data. 



4.5 Parameter estimation 

In the last section of this chapter we make some general remark about the depen- 
dence of the CMB anisotropy spectrum on different parameters and about degen- 
eracies. We start by enumerating the relevant physical processes. 



4.5.1 Physical processes 

• Before recombination, photons and baryons form a tightly coupled fluid which 
performs acoustic oscillations on sub-horizon scales. 

• Depending on the initial conditions, these oscillations are sine waves (isocur- 
vature case) or cosine waves (adiabatic case). 

• After recombination, the photons move along perturbed geodesies, only influ- 
enced by the metric perturbations. 

• Vector perturbations of the metric decay as after creation and their effects 
on CMB anisotropies are negligible for models where initial fluctuations are 
created early, e.g. during an inflationary phase. This is different for models 
which constantly seed fluctuations in the geometry, e.g. topological defects. 

• Tensor perturbations of the metric have constant amplitude on super-horizon 
scales and perform damped oscillations oc once they enter the horizon. 

• Scalar perturbations of the metric are roughly constant if they enter the hori- 
zon only after the time of matter and radiation equality. On scales which enter 
the horizon before equality they are damped by a factor (-2cq/-2in)^, where 2;eq 
and the redshift of equality and of horizon crossing, respectively. 

• Perturbations on small scales, k ^ kx — (i7;,/i/20)(^;dec + l)^-f^o are exponen- 
tially damped by coUisional damping during recombination (Silk damping). 



4.5.2 Scale dependence 

• On large scales (larger than the horizon scale at recombination, £ ^ £h — 
T^I'dui with -Qu = 77doc/x(^o — Vdcc), perturbations are dominated by gravi- 
tational effects: Inflationary models typically lead to k^ — $P(/i;, ?7dcc)) — 
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const, and (H"^) ~ const, on these scales. This imphes a flat "Harrison- 
Zel'dovich" spectrum, 



• On intermediate scales, £h < £ < ^damp ~ 800, CMB anisotropies mainly 
reflect the acoustic oscillations of the photon/baryon plasma prior to recombi- 
nation. The position of the first peak is severely affected by initial conditions 
(adiabatic or isocurvature). For k — 0, the first contraction peak is at about 
i^i^ ~ 220 if the initial conditions are adiabatic, while the first contraction 
peak is at (^^ ~ 350 for isocurvature initial conditions. The amplitude of and 
the distance between the peaks depend strongly on cosmological parameters. 

• On small scales, ^damp < ^, fluctuations are coUisionally damped during recom- 
bination ("Silk damping"). The damping scale depends mainly on flbh and 



4.5.3 The main influence of cosmological parameters 

• Curvature, h'^flf^: 

— Mainly affects the inter-peak distance, Ai, and, for given initial condi- 
tions, the position of the first peak. Positive curvature lowers A£ while 
negative curvature enhances it. 

— Curvature also leads to an integrated Sachs-Wolfe contribution which 
is especially important for k > at very low £. Overall, this leads to 
some enhancement of the Sachs-Wolfe contribution and therefore (after 
normalization to the COBE measurements) to somewhat lower acoustic 
peaks. 

• Bciryon density, pb — Q^bh^ ■ 10~^^g/cm^: 

— A high baryonic density enhances the compression peaks and decreases 
the expansion peaks via the self-gravity of the baryons. 

— It also reduces the damping scale. Ay = l/(adecCT?^e(^dec))) leading to an 
increase in -^damp- 

— Baryons decrease the plasma sound velocity, Cs — V3(l -l- ph/ p-y)~^-i and 
hence prolongs the oscillation period. This increases the spacing between 
acoustic peaks. 

• Cosmological Constant, A = ■ IQ-^^g/cm^: 

The presence of a cosmological constant at fixed Qtot — + delays the 
epoch of equal matter and radiation. During the radiation dominated era, the 




(4.95) 
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gravitational potential is not constant, but decays as soon as a given scale 
enters the horizon. If ?7eq ~ r^dec this induces an integrated Sachs-Wolfe (ISW) 
contribution which boosts mainly the first acoustic peak. also boosts the 
late integrated Sachs-Wolfe contribution. 

• Hubble Parameter, Hq = lOO/i km/(s Mpc): The influence of the Hubble 
parameter is complicated and depends sensitively on the variables which are 
kept fixed during its variation or = h'^Vt,). As one example of its 
influence: for fixed curvature and cosmological constant, lowering the Hubble 
parameter also delays the epoch of equal matter and radiation, r/eq Vdec, 
since 

^eq + 1 = ^ 2.4 ■ lO^Q^h^ . (4.96) 

' 'rad 

Therefore the same type of ISW contribution as for A-models boosts the first 
acoustic peak. 

• Initial conditions: 

— A tensor contribution enhances the large scales fluctuations but not the 
acoustic peaks, thereby lowering their relative amplitude. 

— A "blue" fluctuation spectrum, n > 1, enhances fluctuations on smaller 
scales and raises thereby the acoustic peaks. 



4.5.4 Degeneracy 

One important issue in determining cosmological parameters from CMB anisotropy 
measurements is the choise of good variables, which requires physical insight in how 
anisotropics are influenced. As we have argued before, the Hubble parameter, h 
is not a good variable since its influence is very complicated. It enters the cosmic 
densities p, oc Q^h"^ and the length scales like r/gq or ?7dec- Another limitation for 
parameter estimation from CMB anisotropics is degeneracy. We illustrate here just 
one example. As we have discussed in Chapter 3, the position of the first acoustic 
peak only depends on the sound horizon, = f^'^'"' Cgdrj and the angular diameter 
distance to the last scattering surface, x(^o~''7doc)- The distance between subsequent 
peaks in the CMB power spectrum is proportional to 

^ XijiQ - Vdec) 

In Fig. (left panel) we show lines of constant R = Ai/AiQ in the Qm ^ 
plane. Here AIq = A£o(fiA = = 0) is the value of Ai in a spatially flat universe 
with vanishing cosmological constant. To the right the CMB anisotropy spectra 
for scalar perturbations with fixed index n = 1 and fixed values of the matter 
density and the baryon density Ut- But the cosmological constant and h vary. 
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so that and fim, correspond to the vahies indicated by bullets on the left panel. 
Clearly, for i > 50 these spectra are perfectly degenerate. On the other hand, due 
to cosmic variance, the low £ CMB anisotropics will never be known to very good 
accuracy so that this degeneracy cannot be lifted by CMB anisotropy observations 
alone. Additional data like the supernova type la measurements, observations of 
the galaxy distribution (large scale structure) or CMB polarization are needed. 




Figure 4.4: Left: The lines of constant R are shown in the 0:^-^1^ plane. The 
values f^A, for which the CMB anisotropy spectra are shown right are indicated 
as black dots. Right: Three CMB anisotropy spectra with different values of ^a 
and but fixed R are shown. For £ ~ 50 these spectra are clearly degenerate.The 
solid hne represents a flat model, while the dotted line corresponds to a closed and 
the dashed line to an open universe. 

There are also other degeneracies like the optical depth to reionization and the 
tensor contribution or the scalar spectral index and the tensor contribution. The 
important lesson to learn is that even if the very stringent model assumptions are 
correct, we still need other data to measure cosmological parameters and especially 
we will only feel comfortable with a sufficient amount of redundancy. 
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Chapter 5 

Observations and Results 



In this short, final chapter we want to discuss briefiy the experimental situation 
which is very much in fiow and may have changed considerably already at the mo- 
ment when this review appears. It has been clear for a long time that, if initial 
fiuctuations have led to the formation of large scale structure by gravitational insta- 
bility, they should have induced fiuctuations in the cosmic microwave background 
PB| , Before spring 1992, however, only the dipole anisotropy had been detected 



12, 131. Its value is 




= (1.528 ±0.004) X 10" 



After many upper limits, the DMR experiment aboard the COBE satellite mea- 
sured for the first time convincingly positive anisotropies P]. It found 




- {SOfiKy (5.1) 
on all angular scales 6 > 7°. Many more positive measurements have been performed 



since then. In Fig. |5.1| we just show the COBE DMR results |^ together with the 
three most recent experiments, BOOMERANG MAXIMA- 1 [0 and DASI ^ 

As one sees in this figure, present data, agrees very well with a simple fiat model 
of purely scalar, scale invariant, = 1, adiabatic fiuctuations with cosmological 
parameter uJi, = 0.02, fi^ = 0, Q/^ = 0.7, h = 0.65 which are also preferred from 
other cosmological data. However, the error-bars are still considerable. 

The experiments can be split into three classes: Satellite experiments, balloon- 
borne experiments and ground based experiments. The technical and economical 
advantages of ground based experiments are obvious. Their main problem is atmo- 
spheric fiuctuation. This can be reduced by two methods: 

• Choose a very high altitude and very cold site, e.g. the south pole. Several 
experiments like SP, Python and White Dish have chosen this site. 
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Figure 5.1: The measured temperature anisotropies, i{i+l)Ce are shown in a hn-hn 
plot (left) and in a log-lin plot (right) with the theoretical curve from a standard, 
adiabatic cold dark matter model. The data points shown are those from COBE 
DMR (sohd, magenta, low i), BOOMERANG (solid, red), DASI (dashed, blue) and 
MAXIMA- 1 (dotted, green). 

• Measure anisotropies on small scales, preferably by interferometry (DASI, 
CAT, VSA, Jodrell Bank). 

Balloon-borne experiments flying at about 40km altitude have less problems 
with the Earths atmosphere but they face the following difficulties: 

• They are limited in weight. 

• They cannot be manipulated at will in flight. 

• They have a rather short duration. 

• To secure all the data taken on the balloon, they have to be recovered intact. 

Yet the advantages of overcoming the atmosphere are so significant that many 
groups have chosen this approach, like e.g. MAXIMA-1, TopHat, etc. The BOOMERanG 
experiment combines the two advantages of a cold site and balloon altitude. It has 
performed a long-duration fiight (10 days) on the south pole in December 1998. 

The third possibility are satellite experiments. They avoid atmospheric problems 
altogether, but this solution is very expensive. So far two satellite experiments have 
been launched: COBE in 1989 (NASA mission) and MAP in June 2001 (Microwave 
Anisotropy Probe, a NASA MIDEX mission), one more is planned: PLANCK, an 
ESA medium size mission of the "Horizon 2000" program, to be launched in 2007. 

As I am writing this lines, MAP has safely arrived at its destination, the Lagrange 
point L2 of the sun-earth system. It will perform measurements at five frequencies 
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in the range from 22 to 90 GHz. The instruments of PLANCK will operate at nine 
frequencies, covering 20 to 800 GHz. At low frequencies (below 100 GHz) radio 
receivers are used (so called "HEMTs" , high electron mobility transistors) while the 
high frequency instruments are bolometers. Recent progress in detector technologies 
should enable the two new satellites to perform significantly better than COBE - 
the radio receivers of PLANCK, e.g., are supposed to be 1000 times more sensitive 
than the ones used for COBE, and the angular resolution has improved from seven 
degrees to four arc minutes. For more details also on other experiments see 



http : //astro . estec . esa . nl/PLANCK 



http : / /map . gsf c . nasa . gov 



http : / /www . gsf c ■ nasa . gov/ astro/ cobe/ cobe_home . html 



http : / /spectrum . Ibl . gov/ www/max . html 



http : / /oberon . romal . inf n . it/boomerang/ 



I finish this short chapter with Table ^.1| which shows the ranges for the cos- 
mological parameters Qtot = 1 — ^k, h'^Qt and rig as determined purely by CMB 
anisotropies. Except for the last reference, a purely scalar spectrum of adiabatic 
fluctuations is assumed. The parameter estimation process also assumes 'weak pri- 
ors' on the values of other cosmological observables, like e.g. that the age of the 
Universe be larger than lOGyrs. or 0.4 < h < 0.9. 1 do not comment this table much 
further but refer the reader to the original literature and many improved papers on 
this subject which will appear shortly. 

Clearly, the results shown in Table f).l\ are very consistent. It is interesting to 
note, how the upper limit on the scalar spectral index deteriorates if one allows for 
a tensor component. This is one of the degeneracies in the CMB data which can 
be broken by including large scale structure data in the analysis (see [0). Other 
cosmological parameters are not well constrained by CMB data alone. However, if 
CMB data is combined with SNla and large scale structure data, the error bars 
are significantly reduced and evidence for a non-vanishing cosmological constant 

~ 0.7 becomes very strong (see PO, HO, H2[). 
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Rcf. 


Data 


^tot 




n 


s 


errors 


m 


BOOM and DMR data 


1 02+°-°*^ 

-•^•^^-0.05 


U.UZZ_Q QQ3 


0.96: 


hO.l 
-0.09 


1 


— cr 


errors 




DASI and DMR data 


^■^■-'-0.06 


U.UZZ_Q QQ4 


i.or 


h0.09 
-0.07 


1 


— cr 


errors 




MAX and DMR data 


"•="^-0.16 


U.UOZO„Q QJ25 
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h0.14 
-0.14 


2 


— a 


errors 


il 


all data, no priors 
allows also tensor mode 


1 06+°-^^ 

-•^•^"-0.13 


02+° "^ 


0.93: 


hO.75 
-0.16 


2 


— a 


errors 



Table 5.1: Some results from parameter estimations from recent CMB data alone. 
The errors given are formal 1 or 2-a errors which assume the underlying model to 
be correct and no systematic problems in the data. They are obtained by marginal- 
ization or maximization over all other model parameters. 
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Appendix A 

The Q's from gravitational waves 



We consider metric perturbations which are produced by some isotropic random 
process (for example during inflation). After production, they evolve according to a 
deterministic equation of motion. By reasons of isotropy and due to symmetry, the 
correlation functions of hij(k, rj) have to be of the form 

{hij{k,T])h*i^{k,r]')) = [kikjkikmHi{k,r],r]') + 

{kikdjm + kikmSji + kjki6im + kjkmSii)H2{k, rj, r]') + 
kikj6imH3{k, T], T]') + kikmSijH^{k, r]' , r]) + 
+6ij6imH4{k,r],r]') + {Su^jm + Sim6ji)Hr,{k,r],r]')] .(A.l) 

Here the functions Ha are functions of the modulus = |k| only. Furthermore, all 
of them except H3 are hermitian in t and t'. This is the most general ansatz for an 
isotropic correlation tensor satisfying the required symmetries. To project out the 
tensorial part of this correlation tensor we act on hij it with the tensor projection 
operator. 



T,^ = prnpn _ (l/2)p.^.p-" with P,, = 5,, - k,kj . 

This yields 



(A.2) 




(A.3) 



From Eq. ( [4.17 ), we then obtain 



/AT, AT, .\ 1 f 3 /AT, AT, , A 
(—{n) — in')j = y J d^'x ( — (n,x) — (n',x)J = 




I \ ^ r no rvo 

— 1 / k^dkdil^^ / drj drj' exp^ik ■ n^rjo — r])) exp{—i]<. ■ n{riQ 



v'))- 



(A.4) 
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Here dQ^ denotes the integral over directions in k space. We use the normahzation 
of the Fourier transform 



/(k) 



ci^a;exp(ix- k)/(x) , /(x) 



(2vr) 



(fkexp{-ix ■ k)/(k) 



where V is an (arbitrary) normahzation volume. 

We now introduce the form ( [A.3| ) of < h^'^^h^'^^ >. We further make use of the 
assumption that the perturbations have been created at some early epoch, e.g. dur- 
ing an inflationary phase, after which they evolved deterministically. The function 
H^lk, 7], 7]') is thus a product of the form 



H,ik,ri,ri') = Hik,ri)-H*{k,ri') . 
Introducing this in Eq. ( [A.4| ) yields 
'AT, AT 



(A.5) 



T 



n 



T 

3 



"(nO 



(h) J ^^"^^"^^^ [("^ ■ - 1 + + - W(n • n') + /iV"] ■ 
/ drj drj' H{k,ri)H*{k,ri') exp{ikfi{rio — r])) exp{—ikfi'{rio — r]')) (A. 6) 

•^Vdec -^Vdec 

where /i = (n ■ k) and /i' = {n' ■ k). To proceed, we use the identity 

oo 

exp((zA;/i(r/o - v)) = $^(2r + l)fjr{k{r]o - v))Pr{f^) ■ (A.7) 

r=0 

Here jr denotes the spherical Bessel function of order r and Pr is the Legendre 
polynomial of degree r. 

Furthermore, we replace each factor of /i in Eq. ( |A.6| ) by a derivative of the 
exponential exp{ikfi{riQ — i])) with respect to k{riQ — rj) and correspondingly with /i'. 
We then obtain 



AT, .AT 



T 



n 



T 

3 



n 



(^) X](2r + l)(2r' + l)i(^-'-') J edkdn^^PMPr'ifi') 
2(n ■ n )^ / dr]dr]'jr{k{r]o - r]))jr'{k{r]o - r]'))H{k,r])H*{k,r]') 



dr]dr]'[jr{k{r]o - r]))jr'{k{r]o - r]')) +£{k{r]o - v))jr'{k{r]o - r]')) + 
Ukivo - vM'iHVo - V')) - J'mvo - vM'iHVo - V'm{k, r^)H*{k, V) 



4(n-n') / dr]dri'fMr]o-v))j'Akim-v'))Hik,r])H*{k,ri') 
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Here only the Legendre polynomials, Pr{fi) and Pr'ifJ^') depend on the direction 
k. To perform the integration dQ^., we use the addition theorem for the spherical 



harmonics Yr. 



An 



The orthogonality of the spherical harmonics then yields 

(2r + l)(2r' + l) ^ dn^Pr{ii)Prifi') 

r 

s=~r 

4:7r6rr'Pr(n ' h') . 



(A.9) 



(A.IO) 



In Eq. ( |A.8| ) the integration over dQ^. then leads to terms of the form (n-n')Pr(n-n') 
and (n ■ n')^P,.(n ■ n'). To reduce them, we use 



r + 1 r 

xPr(x) = Pr+l H Pr-l 

Applying this and its iteration for x'^Pr{x), we obtain 



(A.ll) 



AT AT* 



^^(2r + l) J k^dk j dridri'H{k,ri)H*{k,ri')^^ 
2(r + l)(r + 2; 



-Pr 



(2r + l)(2r + 3)^''^^ ' (2r - l)(2r + 3) ^ ' ' (2r-l)(2r + l) 
jr{k{r]o - r]))jr{k{r]o - r]')) - Pr[jr{k{r]o - v)jr{k{Vo - v')) 

+jr{k{V0 - V'))jriKVO - V)) - JriKVO " V))jr'iKVO " V'))] 



-Pr 



2r(r 



Pr 



r-2 



X 



-4 



r + 1 
2r + 1 



Pr+1 + 



2r + 1 



Pr^l 



(A.12) 



where the argument of the Legendre polynomials, n- n', has been suppressed. Using 
the relations 

for Bessel functions, and its iteration for j" , we can rewrite Eq. ( |A.12| ) in terms of 
the Bessel functions jr-2 to jr+2- 

We now insert the definition of Cg: 



/AT, AT , 
.n ■ — — n 



\ T 



T 



(n-n')=cos ( 



1 

An 



S^(2£+ l)QP^(cos( 



(A. 14) 
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and compare the coefficients in Eqs. ( [A.12| ) and ( |A.14|) . We obtain the somewhat 
lengthy expression 



2 



dkk^ / dr]dr]'H{k,ri)H*{k,ri')\ji{k{rio-ri))ji{kirio-ri')) X 



1 ^2(2£2 + 2£-l) ^ (2£2 + 2£-l)2 



(2£-l)(2£ + 3) (2£-l)(2^ + 3) (2£- l)2(2£ + 3)2 
4£3 4(£+l)3 



(2£- 1)2(2£+ 1) (2£+ l)(2£ + 3)2^ 

[j£(^(^o - v))je+2{Hvo - v')) + JwiHvo - v))jfXHVo - v'))] X 
1 ,^2(£ + 2)(£+l)(2£2 + 2£-l) ^ 2(£ +!)(£ + 2) 8(£ + l)2(£ + 2) 



2/ + 1V (2£- l)(2£ + 3)2 (2£ + 3) (2£ + 3)2 

'2£(£-l)(2£2 + 2£-l) 2i{£~l) 8£2(£-l) 



2/ + 1V (2£- l)2(2£ + 3) (2£-l)(2 (2£ - 1)2 
2(£ + 2)(£ + l) 4(£ + l)(£ + 2)2 ^ (£+i)2(£ + 2)2 



(2£+l)(2£ + 3) (2£+ l)(2£ + 3)2 (2£+ l)2(2£ + 3)2 
+je-2{k{rio - v))ji-2{k{r]o - V)) x 



(2£-l)(2£ + l) (2£- 1)2(2£+ 1) (2£-l)2(2£ + l)^ 

An analysis of the coefficient of each term reveals that the curly bracket in this 
expression is just 

and the result is equivalent to 

Ci = ^ j dkk^\I{i,k)\H{i-l){i +!){£ + 2) , (A.16) 

with 

/(t^^ /'%,gfat y;y<''°- ''» . (A.17) 
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Appendix B 

Boltzmann equation and 
polarization 



In this appendix we derive the Boltzmann equation taking into account polarization, 
and we write it as a hierarchy of equations using an orthonormal expansion in the 
space of photon directions. Up to the collision term, the Eqs. (^4.65|) , (^4.66|) and 



( [4.67| ) are still valid. We first re-derive the collision term taking into account the 
polarization dependence of Thomson scattering. 

Just before the process of recombination during which the fluid description of 
radiation breaks down, the temperature is ~ 0.4 eV. The electrons and nuclei 
are non-relativistic and the dominant collision process is non-relativistic Thomson 
scattering. 

Thomson scattering depends on the polarization of the incoming radiation field. 
We describe the polarization state of the radiation field by the Stokes parameters |T^ , 



20, 21, IS 



For a harmonic electro-magnetic wave with electric field 

E(x, t) = (ei^i + €2^2) e*^^ ■ , (B.l) 

where n, ei and €2 form an orthonormal basis and the complex field amplitudes are 
parameterized as Ej = ajc''^^ , the Stokes parameters are given by 

I = aj + aj (B.2) 

Q = aj-al (B.3) 

U = 2aia2Cos{S2- 5i) (B.4) 

V = 2aia2sin((52 (B.5) 

/ is the intensity of the wave (whose perturbation Ai has interested us so far), while 
Q is a measure of the strength of linear polarization (the ratio of the principal axis 
of the polarization ellipse). U and V give phase information (the orientation of 
the ellipse). For non-relativistic Thomson scattering V is completely decoupled and 
(since it vanishes at early times) is therefore never generated. 
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As Q and U vanish in the background, perturbations cannot couple to them 
(since such terms are 2nd order), and the Liouville equations for Q and U become 
(neglecting scattering and spatial curvature) 

d^r,iQ;U)+in'h{Q;U) = 0. (B.6) 

The differential cross section of Thomson scattering for a photon with incident 
polarization e(j) scattering into the outgoing polarization e^s) = e' is fl^ 



n 




Figure B.l: Definition of the angles and vectors for Thomson scattering in the 
(n, 62) plane. 



It is often convenient to introduce the two 'partial' intensities Ji 



{I — Q)/2. A wave scattered in the (n, £2) plane (see figure |B.1| ) by 



and I2 = a 
an angle 6 has the intensities 



is) 



3aT (i) 
"-'1 



r(«) 



»7r 

^/fcos^e. 



or, expressed in terms of the Stokes parameters, 

3(Tt / 1 + cos2 e sin^ 6 



lis) 



167r 



sin^ 6 1 + cos^ 



(B.8) 



:b.9) 



A rotation in the (ei, 62) plane doesn't change the intensity of the wave, but it 
changes Q and U to 



Q' = Qcos(20) + f/sin(2(/)) 
U' = -f/sin(20) +gcos(20) . 



(B.IO) 
(B.ll) 



To determine the cross section that a given 'initial' wave 
{I^^\Q^^\U^'^^) propagating in direction n be scattered into a wave {I^'^\Q^^\U^'^^) 
with direction n', we need to go through the following steps (we will use the plane 
{y, z) as reference plane, see figure ( p.2| ) for definitions of the angles and vectors): 
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Figure B.2: Definition of tlie angles and vectors for Thomson scattering in tfie 
general case. The polarization vectors are oriented like in figure |B.l. 



Rotate around n such that the plane (n, n') turns into the plane (n^;). One 
needs to apply the rotation ( |B.1(J| , [B!TT| ) for = a to the Stokes parameters. 

Rotate the new plane (n, n') around z into the reference plane {y,z). This 
operation does not influence the incoming Stokes parameters.. 



3. Now we are in the known case of ( p.8|) and (|B.9|) . Hence we can apply the 
scattering matrix. 

4. We then rotate the scattering plane back around z into the old (z, n') plane. 
This does not change the scattered Stokes parameters. 

5. Finally we rotate around n' by the angle a' to reach the original state. To 
do this, we have to apply the rotation matrix ( p.lCI| , |RTTD again, but now for 
= a'. 

Following the steps outlined above, we recover the scattering matrix in the basis 
{Ii,l2,U) given in equations ( |B.13|) - (|B.16| ) (see also [1^). V is completely de- 
coupled from the other parameters and follows an evolution which is independent 
of the rest. Hence by starting with V{t -C tdec) = it will stay zero and can be 
neglected. The scattering matrix P, which determines the (non vanishing) scattered 
Stokes parameters from the initial ones, 

if I (B.12) 




47r 



'2 
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is then given by 



P 



(B.13) 



where 



p(0) 



p{i) 




4/XyU,' cos(0' — 


— 4/i' sin(0' — 



2/isin(0'- 


2cos('(/)'- 




p(2) ^ :d 

4 



/x2,i'2cos[2(0' - (f>)] 
-fj.'^ cos[2((/.' - (/))] 
-2^t/i'2 sin[2(0' - </))] 



-fj.^ cos[2(<j>' - (/>)] ^j.^ ii' sm[2{<f>' - 4>)] 
cos[2{</)' - (/))] -fi' sin[2(fli' - 0)] 
2/i sin[2((/)' - </))] 2/i/i' cos[2((/.' - </>)] 



(B.14) 

(B.15) 
(B.16) 



As we are in an isotropic situation, we will perform all the calculations in a 
special coordinate system with k || i and n, n' as in Fig. |B.2| . Clearly the results are 
independent of this coordinate choice. 

The matrix R connecting (/i, I2, U) with (/, Q, U) is given by 



1/2(7 + Q) 
1/2(7 -Q) 
U 





:b.i7) 



To calculate the collision term including polarization , we change into the (7i, I2) 
basis . For each of the two intensities A G {1,2} we then have the collision term 
given by 



C[f 



{A)l 



(A) 



(B.18) 



dr] dr] 

where and denote the distribution of photons in the polarization state A 
scattered into respectively out of the beam due to Compton scattering. 

In the matter (baryon/electron) rest frame, which we indicate by a prime, we 
know that 

^(p,n) = f^^yip',n')Ps\n,n')dn' , 

where He denotes the electron number density and is the 2x2 upper left corner of 
the normalized Thomson scattering matrix ( p.l3|) . In the baryon rest frame which 
moves with four velocity u, the photon energy is given by 

p' = p^uf" . 

We denote the photon energy with respect to a tetrad adapted to the slicing of 
space-time into {77 = constant} hyper-surfaces by p : 



p = p^n"^ , 



with n = a-^[{l-A)d^ + B'di 
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The lapse function and the shift vector of the shcing are given hj a = a{l + A) and 
(3 = —B^di . In first order, 

Pq = ap{l + A) — apriiB^ , 

and to zeroth order pi = aprii. Furthermore, the baryon four velocity has the form 
= a~^{l — A) , = up to first order. This yields 

p' = p^u^ =p{l + n,{v' -B')) . 

Using this identity and performing the integration over photon energies, we obtain 



drj' 



Ait 



Photons which are scattered leave the beam, with the probability given by the 
Thomson cross section (see e.g. 



df. 



(A) 



dt' 



so that we finally have 



dp 



df^ df_ 
dt' 



(A)' 



dt' 



p 



'-aTne[Ani{v' - B') 



,(A) 



(B.19) 



By setting C^^'^ = + C(2) and = - C(2) we transform the collision 
integral back to the normal stokes parameters. The term for U has the same form 
as the one for Q, just with the corresponding matrix elements in the integral. The 
Boltzmann equation then finally becomes (setting S = {Ai,Q,U) and for the flat 
case, K = 0): 



M + ifikM = Aipk{^ - ^ + n™S(;^)) + An^rf'H, 



ml 



+aaTne 



Q + ifikQ = aaxrif. 



II + ipkU = aarUe 



-Q + I dn'p^s'^ 



-u + I dn'p^s'^ 



(B.20) 

(B.21) 
(B.22) 
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where we have to use the scattering matrix transformed into the {A4, Q, U) basis, 



P = Ps + Pv + Pt with (B.23) 
Ps = R'^P^^^R (B.24) 

{i-/.2)(i_3^,2) 3(i-m2){i-m'^) (B.25) 







V 



1 - fi^^l- ^I'^R-^P^^^R (B.26) 

fifi'C fifi'C -fiS (B.27) 
fj.'S fx'S C I 

Pj, = R-^P^^^R (B.28) 

3/ {1-|.2)(1-A.'2)Ct -(l~tJ?){l + tl'^)CT 2(1-|.2V5t \ 
= - -(1 + |.2)(1-m'2)Ct {1+a.2)(i + ^'2)Ct -2(1 + m2)m'St (B.29) 

8\ -2^1(1 - At'2)5T 2/i(l +/i'2)ST WCt / 

with C = cos(0 — 0')) "S" = sin(0 — 0') and 

Ct = cos(2(0 - 0')), St = sin(2(0 - 0'))- The parts Ps, Py, Pt of P describe the 
scattering of the scalar, vector and tensor contribution to S respectively. 

The functions Ai, Q and U depend on the wave vector k, the photon direction 
n and conformal time 7]. We choose for each k-mode a reference system with z-axis 
parallel to k. For scalar perturbations we achieve in this way azimuthal symmetry 
— the right-hand side of the Boltzmann equation and therefore also the brightness 
Ad^^"^ depend only on /i = (k-n) and can be expanded in Legendre polynomials. The 
right-hand side of the Boltzmann equation also determines the azimuthal dependence 
of vector and tensor perturbations. One can continue with this approach, but the 
resulting equations for Q and U and especially also their power spectra depend 
explicitly on the coordinate system. Therefore, we prefer an approach which is 
inherently covariant under rotation. 



B.l Electric and magnetic polarization 



Since the Stokes parameters explicitly depend on the coordinate system, and Eqs. ( |B.21| ) 



and ( |B.22| ) transform rather complicated under rotations of the coordinate system. 
A method to characterize CMB polarization due to non- relativist ic Thomson scat- 
tering which is more convenient than the Stokes parameters since its transformation 

A 
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properties are very simple has been developed a couple years ago [0, |49| , 
detailed derivation of this method goes beyond the scope of this report. Here we 
just repeat the definitions and the main results. We set 

r = {M,Q + tU,Q ~iU) (B.30) 

In terms of this vector the collision integral above can we written (in vector form) 
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as 



c\r] 



[-T+ (^i- j dfi'A<' + (n-v,),0,0^ 
^ E / dfi'P(™)(n,n')T'] (B.31^ 

m=-2 



From Eqs. ( [B.13| ) to ( [B.17| ) one can determine the scattering matrix for the vector 



p{m) 



/ 3 \rml-V'r. 



62^2 2-^2 
6 21 2 -21 2 



0^21 2 2^2 



3 -an™'- 



2I2 



\ 



J 



^B.32) 



where sY^^' = sYl^*{n') and gY^^ are the spin-weighted spherical harmonics ^ 



We now decompose the Fourier components of the temperature anisotropy M. 
and the polarization variables E and B as 



i m=-2 



[B.33) 

;b.34) 



m=-2 



Here m = is the scalar mode, m = ±1 are the vector and m = ±2 are the tensor 
modes. The functions are closely related to the spin weighted harmonics sYf^: 



sGTin) 



4tt 



.Yr(n) 



2i+l' 

The evolution equations in term of these variables can be given in the following 



compact form 



Mi^^ - k 



(B.35) 



E 



(m) 



2£ 



I 771V'"; 



(»n) 2m ^(m) 2l^T+l jp{m) 



+ 11 



2£ + 3^^+^ 



-Up 



B^r^ - k 



21-1 
—rieaTaB) 



2 (B.36) 

(m) _ 2^^+i j^(m) 



+ 1 



2£ + 3" 



'i+i 



(B.37) 
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where we have set 

AO) _ ^ ^_^^(o) c(i) 



5f = ne(TTaC(o), S\^> = ne(7Ta4o;6, (B.38) 

and C("^) = 15 [^^2""^ - v^^^:^"'^]. The couphng coefficients are 



Note that for scalar perturbations, m = 0, i?-polarization is not sourced and we 
have Bf = 0. 

Finally we want to connect the intensities Aif^^ with the more familiar expansion 
of the scalar (S), vector (V) and tensor (T) contributions to the brightness function 
in terms of Legendre polynomials. Usually one sets 

Here Ad'^^^ only depends on = (n ■ k)/A; and the n-dependence of M.^^^ and M.^'^^ 
can be written as 

M^^\fi,(p) = x/l-Ai' \M^'\fi) cos (^+M^Pifi) sin cp\ (B.39) 
M^^\li, 0) = (1 - n^) \m^+'^ cos(20) + ) sin(20)] , (B.40) 
where (j) is the azimuthal angle in the plane normal to k. Each of the functions 

( S V 

is now expanded in Legendre polynomials 

t 

The coefficients cr^'^'^'^^ are then related to Al^"*^ via the identities 

Mf = (2£+l)af^ (B.42) 



(±2) _ _ (£ + 2)! 1 I (^) 

~ y (£-2)!^2£ + 3 (2£-l)(2£ + 3) Ti.^ 

+^'-n.l2] , (B.44) 

We do not repeat this correspondence for the Stokes parameters Q and U since 
it is rather complicated and not very useful as it depends on the coordinate system 
chosen. 



where 
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B.2 Power spectra 

In the previous appendix and in Chapter 4 we have derived the expression for the 
CMB anisotropy power spectrum for scalar and tensor perturbations. Here we give 
the general expression for scalar, vector and tensor fluctuations, polarizations and 
cross correlations. To make contact with the results derived before, one has to 
use Eqs. ( [B.41| , P^^ ) and (|B.44[ ) and neglect the collision term in the Boltzmann 



equation. 

We expand the present CMB anisotropics and polarization in spherical harmon- 
ics: AT(n, rjo) /Tq = J2em '^^m^'"(n). The coefficients aim are random variables with 
zero mean and rotationally invariant variances, Cg = {\ aim P)- The correlation 
function of the anisotropy pattern then has the standard expression: 



ni)^(n2)\ = ^J2(^i+l)C,Pi{cose) (B.45) 



ST 

where cos^ = rii ■ n2 and (■ • ■) denotes ensemble average. We use the Fourier 
transform normalization 

/(k) = /(x) exp(zk ■ ^)d'x , (B.46) 

with some normalization volume V. Using statistical homogeneity we have 

5T. .5T. \ 1 f 3 /5T. .5T. 
— (ni)7;r(n2) ) = t7 d ^ \ 7;r(x, ni) — (x, ns, 



To' ''To' V J \To' ' "To 

1 f /6T „ JT. 



(2vr; 



J rf3fc^^(k,ni)^(k,n2)^ . (B.47) 



Inserting our ansatz ( |B.41| ) for ^ = jAi^ and using the addition theorem for spher- 



ical harmonics 



Pt{ni .112) = ^ E^YemMYimM, we find 



^(n,)^(n,)) = i- Y: i-lf'''Y,Un.)Y;Un2) 



To ' 'To' ' Stt 

1,1' ,m,m' 



X j edkdQ^Y;^{k)Yem'{k){a,a*,){k) 
= ^J2^2e +l)Pe{ni- 112) f k^dk{aea'^){k) , (B.48) 

from which we conclude 

Ct''''^'^ = ^ fedk{\af\to,k)\'), (B.49) 
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where the superscript indicates that Eq. ( |B.49| ) gives the contribution from scalar 
perturbations and -^-^ means that it is the contribution to the intensity perturba- 
tion. 

The QQ, UU, A4Q, M.U and QU correlators depend with the Stokes parameters 
on the particular coordinate system chosen. It is much more convenient to express 
the polarization power spectra in terms of the variables E and B which are indepen- 
dent of the coordinate system. Furthermore, since B is parity odd, its correlators 
with M. and E vanishes. One finds the simple general expression [^| 

(2£ + lfCf^^'^'> = / k^dkXf'^Y^'^'^* , (B.50) 

where rim = 1 for rn = Q and ra^ = 2 for m = 1,2, accounting for the number of 
modes. 
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